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I 


INTRODUCTION 


I.  INTRODUCTION 


A mathematical  model  suitable  for  simulation  of  bridge  collision  inci- 
dents involving  inland  waterway  tows  "has  been  developed,  as  described  in 
the  accompanying  report^.  The  mathematical  model  has  been  implemented  in 
a FORTRAN  computer  program,  and  is  designed  to  run  in  a time  sharing  or 
dedicated  processer  environment,  with  interaction  from  the  user  throughout 
the  simulation.  The  input  data  file  format  is  simple  and  consistent,  and 
all  interaction  with  the  program  is  in  response  to  on-line  prompts.  The 
structure  of  the  input  data  file  is  described  in  the  next  section,  and  a 
discussion  of  how  to  run  the  program  and  interpret  the  output  is  given  in 
the  following  section.  The  discussion  is  presented  in  relation  to  a sample 
case,  and  the  corresponding  input  and  output  formats  are  illustrated  in 
the  Appendices.  A program  listing  is  also  appended  to  this  report. 


II.  INPUT  DATA  FILE  DESCRIPTION 

The  input  data  file  format  has  been  structured  to  be  as  convenient  as 
possible  from  the  users  prospective.  Each  line  of  the  data  file  is  coded, 
to  identify  the  type  of  data  contained  on  that  line,  and  all  data  lines  are 
arranged  in  the  same  format.  Therefore,  the  user  can  check  very  easily 
whether  or  not  all  of  the  required  data  elements  have  been  included  in  the 
data  file. 

The  program  is  structured  to  enable  several  simulation  cases  to  be  ex- 
ecuted consecutively  from  the  same  data  file.  To  execute  several  cases 
consecutively  it  is  necessary  only  to  specify  the  data  values  for  each  case 
which  differ  from  the  preceeding  case.  Furthermore,  the  requirements  for 
ordering  the  lines  in  the  data  file  are  minimal;  the  first  line  of  the 
file  must  contain  the  text  of  a heading  that  appears  on  the  output,  the 
last  line  of  data  for  each  case  must  be  a case  delimiter,  and  the  last  line 
of  the  data  file  must  be  a run  delimiter.  The  lines  of  data  within  each 
case  may  be  arranged  in  any  order.  Cases  are  executed  sequentially,  in  the 
order  in  which  they  appear  in  the  data  file. 
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The  program  assumes  the  first  line  in  the  data  file  has  the  format  (20A4). 
All  remaining  lines  are  assumed  to  have  the  format  (A4,  IX,  11,  13,  IX,  7F10.5). 
The  characters  and  values  in  each  field  are  read,  one  line  at  a time,  and 
assigned  to  the  following  temporary  variables: 


Column 

Type 

Name 

1-  4 

Alpha 

LABEL 

6 

Integer 

ITAG 

7-  9 

Integer 

JTAG 

11-20 

Fixed 

TEMP (1) 

21-30 

Fixed 

TEMP(2) 

31-40 

Fixed 

TEMP(3) 

41-50 

Fixed 

TEMP (4) 

51-60 

Fixed 

TEMP (5) 

61-70 

Fixed 

TEMP (6) 

71-80 

Fixed 

TEMP(7) 

The  value  of  LABEL  in  each  line  specifies  the  type  of  data  contained  in 


that  line.  The  permissible  character  strings  for  LABEL,  and  the  type  of 

data  to  which  each  corresponds,  are: 

BARG  - data  describes  the  characteristics  of  the  barge  array 
or  the  tugboat, 

INIT  - data  describes  the  initial  conditions  (velocity  and 
orientation)  of  the  flotilla, 

ROUT  - data  describes  a segment  of  a route, 

COEF  - data  specifies  values  for  hydrodynamic  coefficients, 

STER  - data  specifies  steering  criteria  parameters, 

TRIP  - case  delimiter;  causes  execution  of  a case  to  begin, 

QUIT  - run  delimeter;  causes  program  execution  to  terminate. 

A typical  data  file  may  consist  of  a heading  (TITLE)  line,  five  BARG 
lines,  two  INIT  lines,  a STER  line,  four  COEF  lines  and  a number  of  ROUT 
lines.  The  input  for  the  first  case  is  terminated  with  a TRIP  line.  Sub- 
sequent cases  need  contain  only  the  lines  that  are  changed  from  the  preceed- 
ing  case,  and  are  each  (including  the  last)  terminated  by  a TRIP  line.  The 
last  line  in  the  data  file  is  a QUIT  line.  A typical  input  data  file  is 
shown  in  Appendix  A.  The  content  of  each  type  of  data  line  is  discussed 
below. 
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Heading 

The  first  line  of  the  data  file;  its  content  is  printed  as  a heading 
on  the  output. 


Barge  and  Tugboat  Characteristics 

Five  lines  of  data  are  required  to  describe  the  dimensions  and  physical 
characteristics  of  the  flotilla,  its  rudder  and  its  propeller. 


Name  Units  Entry 


LABEL  BARG 

ITAG  1 

JTAG  NBARW,  number  of  barges  wide 

TEMP(l)  ft  TOWWID,  width  of  barge  array 

TEMP(2)  ft  BOTWID,  width  of  tugboat 

Name  Units  Entry 


LABEL  BARG 

ITAG  2 

JXAG  NBARL,  number  of  barges  long 

TEMP (1)  ft  TOWLEN,  length  of  barge  array 

TEMP (2)  ft  BOTLEN,  length  of  tugboat 

TEMP (3)  ft  CGTOW,  LCG  of  barge  array 

TEMP (4)  ft  CGBOT , LCG  of  tugboat 

TEMP (5)  ft  TOWK,  gyradius  of  barge  array 

TEMP (6)  ft  BOTK,  gyradius  of  tugboat 

Name  Units  Entry 


LABEL  BARG 

ITAG  3 

JTAG  0 

TEMP(l)  ft  TOWDRF,  draft  of  barge  array 

TEMP (2)  ft  BOTDRF,  draft  of  tugboat 

TEMP (3)  TOWBC,  block  coef.  of  barge  array 

TEMP(4)  BOTBC,  block  coef.  of  tugboat 

Name  Units  Entry 


LABEL  BARG 

ITAG  4 

JTAG  NPROP , number  of  propellers.  If 


NPROP > 0 four  quadrant  operation 
is  assumed.  If  NPROP  <0,  NPROP  = 
| NPROP | and  only  first  quadrant 
operation  is  possible 


TEMP(l) 

hp 

SHP,  horsepower  per  shaft 

TEMP (2) 

rpm 

RPMAX,  maximum  propeller  speed 

TEMP(3) 

ft 

OFSET , shaft  offset  from  centerline 

TEMP (4) 

ft2 

ARUD,  rudder  area 

TEMP (3) 

rad 

DELMAX,  maximum  rudder  angle 

TEMP (6) 

rad/sec 

DLDTMX,  maximum  rudder  rate 
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Name  Units  Entry 


LABEL  * BARG 

ITAG  5 

JTAG  NPROP 

TEMP(l)  ft  DPROP,  propeller  diameter 

TEMP(2)  ft  PITCH,  propeller  pitch 

TEMP(3)  ARAT,  propeller  blade  area  ratio 

TEMP (4)  WFRAC,  wake  fraction 

TEMP (5)  TDDUC,  thrust  deduction  factor 


Initial  Conditions 

One  line  is  required  to  specify  the  initial  velocity  and  orientation 
of  the  flotilla  as  it  enters  a curved  segment.  Another  line  specifies  the 
integration  step  size  and  error  control  parameters. 


Name 

Units 

Entry 

LABEL 

— 

INIT 

ITAG 

— 

0 

JTAG 

— 

INSEG,  initial  segment  number 

TEMP(2) 

ft/sec 

SPDIN,  initial  speed  of  flotilla 

TEMP (4) 

rad 

GAMIN,  initial  velocity  vector  angle 

TEMP (5) 

rad 

HEADIN,  initial  yaw  angle 

TEMP  (6) 

rad /sec 

CDOTIN,  initial  yaw  rate 

TEMP  (7) 

DRADIN,  initial  radial  offset  factor 

Name 

Units 

Entry 

LABEL 

INIT 

ITAG 

1 

JTAG 

NCUTS,  integration  step  size  limiter 

TEMP (1) 

sec 

STEP,  integration  interval 

TEMP (2) 

sec 

FIRSTP,  nominal  integration  step  size 

TEMP(3) 

EPS,  relative  error  limit 

TEMP  (4) 

AB,  absolute  error  limit 

Route  Description 

One  line  contains 

the  number  of 

segments,  NSEG,  used  to  describe  a 

route.  An  additional 

number,  NSEG, 

of  lines  contain  the  characteristics 

of  each  segment. 

Name 

LABEL 

ITAG 

JTAG 


Units 


Entry 

ROUT 

0 

NSEG,  number  of  last  route  segment 
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Name 

Units 

Entry 

LABEL 

ROUT 

ITAG 

ISUB,  sub-segment  number 

JTAG 

ISEG,  segment  number 

TEMP  (1) 

deg 

SANG  (ISUB,  ISEG),  angular  extent  of 
subsegment  (ISUB  - 1) , if  ITAG  > 1 

mph 

SWIND  (ISEG),  wind  speed,  if  ITAG  = 1 

TEMP (2) 

ft 

SDAT  (ISUB,  ISEG,  1),  radial  distance 
to  inside  wall  of  channel 

TEMP(3) 

ft 

SDAT  (ISUB,  ISEG,  2),  radial  distance 
to  outside  wall  of  channel 

TEMP (4) 

ft/sec 

CUR  (ISUB,  ISEG,  1),  current  speed 
at  inside  wall  of  channel 

TEMP(5) 

ft/sec 

CUR  (ISUB,  ISEG,  2),  current  speed 
at  outside  wall  of  channel 

TEMP (6) 

ft/sec 

CUR  (ISUB,  ISEG,  3),  cross  current 
speed 

TEMP (7) 

deg 

DWIND  (ISEG),  wind  direction  if 

ITAG  = 1 

Hydrodynamic  Coefficients 

As  many  as  four 

lines  are  required  to  specify  values  for  the  hydrodynamic 

coefficients.  If  a 

bowthruster  is 

included  in  the  flotilla,  three  additional 

lines  are  required 

to  satisfy  its 

characteristics . 

Name 

Units 

Entry 

LABEL 

COEF 

ITAG 

1 

JTAG 

1 

TEMP(l) 

— 

AA(2),  coefficient  a2 

TEMP (2) 

AA(3),  coefficient  a3 

TEMPO) 

AA(9),  coefficient  ag 

TEMP (4) 

AA(10),  coefficient  a^g 

Name 

Units 

Entry 

LABEL 

COEF 

ITAG 

1 

JTAG 

2 

TEMP(l) 

AA(1),  coefficient  a^ 

TEMP (2) 

AA(5),  coefficient  a^ 

TEMPO) 

— 

AA(6),  coefficient  af, 

TEMP (4) 

AA(7),  coefficient  a -j 

Name 

Units 

Entry 

LABEL 

COEF 

ITAG 

2 

JTAG 

— 

1 

TEMP (1) 

BB(2),  coefficient  b_ 

TEMP (2) 

BB(3),  coefficient  b3 

TEMPO) 

BB(9),  coefficient  bg 

TEMP (4) 

BB(10),  coefficient  b^g 

L 
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Name 

LABEL 
ITAG 
JTAG 
TEMP  (1) 
TEMP (2) 
TEMP(3) 
TEMP (4) 

Name 


Units 


LABEL 

ITAG 

JTAG 


Units 


coefficient 
coefficient  b5 
coefficient  b^ 
coefficient  b., 


Entry 

COEF 
2 
2 

BB(1), 

BB(5) , 

BB  (6) , 

BB(7) , 

/ 

Entry 

COEF 

4 

NBTSPD;  the  number  of  data  points  des- 
cribing the  bowthruster  performance 
curve;  if  NBTSPD  = 0,  no  bowthruster 
is  included;  required  only  if  one 
was  included  in  a previous  case. 


TEMP(l) 

lb 

BTMAX,  maximum  bowthruster  force 

TEMP (2) 

ft 

BTPOS,  position  of  bowthruster 

TEMPO) 

— 

BTGAIN,  lateral  thrust  command 

P 

Name 

Units 

Entry 

LABEL 

COEF 

ITAG 

5 

JTAG 

1 

TEMP(l) 

ft/sec 

BTSPD(l),  speed  for  data  point 

1 

TEMP  (2) 

ft/sec 

BTSPD(2),  speed  for  data  point 

2 

TEMPO) 

ft/sec 

BTSPDO) , speed  for  data  point 

3 

TEMP (4) 

ft/sec 

BTSPD(4),  speed  for  data  point 

4 

TEMP (5) 

ft/sec 

BTSPD (5),  speed  for  data  point 

5 

TEMP (6) 

ft/sec 

BTSPD(6),  speed  for  data  point 

6 

TEMP(7) 

ft/sec 

BTSPD(7),  speed  for  data  point 

7 

Name 

Units 

Entry 

LABEL 

COEF 

ITAG 

5 

JTAG 

2 

TEMP (1) 

lb 

BTHRUS (1) , thrust  at  BTSPD  (1) 

TEMP (2) 

lb 

BTHRUS (2) , thrust  at  BTSPD  (2) 

TEMPO) 

lb 

BTHRUS (3) , thrust  at  BTSPD  (3) 

TEMP (4) 

lb 

BTHRUS (4) , thrust  at  BTSPD  (4) 

TEMP (5) 

lb 

BTHRUS (5) , thrust  at  BTSPD  (5) 

TEMP (6) 

lb 

BTHRUS (6) , thrust  at  BTSPD  (6) 

TEMP (7) 

lb 

BTHRUS (7) , thrust  at  BTSPD  (7) 

Steering  Criteria 

One  line  is  required  to  specify  values  for  the  rudder  control  parameters. 
This  data  is  not  required  for  solution  by  interpolation. 
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Name  Units  Entry 


LABEL  STER 

ITAG  • 0 

JTAG  0 

TEMP(l)  — BOWCLR,  initial  bow  clearance  factor 


Case  Delimiter 

One  line  is  required  to  initiate  execution  of  a case  and  to  specify  the 
solution  option  selected.  This  line  must  immediately  follow  the  last  line 
of  data  for  each  case. 

Name  Units  Entry 

LABEL  TRIP 

Run  Delimiter 

The  last  line  of  the  data  file;  it  causes  the  program  execution  to 


terminate. 

Name  Units  Entry 

LABEL  QUIT 


III.  PROGRAM  EXECUTION  AND  OUTPUT  DESCRIPTION 

Prior  to  execution  of  the  computer  simulation,  it  is  necessary  to  set  up 
an  input  data  file,  which  will  specify  the  physical  and  hydrodynamic  charac- 
teristics of  the  flotilla,  its  initial  position,  orientation  and  velocity, 
and  the  physical  characteristics  of  the  channel.  The  format  of  the  input 
data  file,  and  the  ordering  of  the  variables  within  each  line,  is  described 
in  the  preceeding  section.  The  data  requirements  in  most  cases  are  entirely 
self  explanatory,  however,  certain  conventions  and  options  which  require  sup- 
plementary explanation  are  described  in  greater  detail  below. 

In  the  BARG  5 data  line,  the  second  integer  entry,  JTAG,  specified  the 
number  of  blades  per  propeller.  However,  if  this  entry  is  made  negative, 
the  absolute  value  will  be  assumed,  but  it  will  serve  as  a flag  to  direct 
the  program  to  use  an  alternate  data  file  for  the  propeller  characteristics. 

II 

I ' 
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The  normal  data  file  has  four  quadrant  data  for  a specific  propeller  configura- 
tion (blade  area  ratio,  pitch/diameter  ratio,  etc.)*  The  alternate  data  file 
has  a more  generalized  structure,  whdre  blade  area  and  pitch/diameter  can 
assume  any  values,  however,  the  data  is  only  applicable  to  ahead  operation 
with  forward  propeller  rotation.  The  option  is  desirable,  however,  if  the 
effect  of  variations  in  propeller  characteristics  are  to  be  evaluated. 

The  physical  characteristics  of  the  channel  are  described  by  the  set  of 
ROUT  data,  as  described  previously.  The  sign  convention  assumed  for  the  cur- 
rent direction  is  that  the  current  runs  opposed  to  the  barge  direction;  that 
is,  for  a segment  where  aQ  > 0,  the  current  runs  clockwise  around  the  origin 
of  the  segment.  When  aQ  < 0,  the  current  runs  conterclockwise  around  the 
origin  of  the  segment.  The  cross  current  is  assumed  to  run  radially  outward 
away  from  the  origin  of  the  segment.  To  reverse  any  of  these  directions, 
simply  specify  the  appropriate  velocity  component  to  be  negative.  Note  par- 
ticularly that  the  current  along  the  inside  and  outside  wall  of  the  channel 
may  have  opposed  signs,  thereby  representing  a current  reversal  or  eddy. 

If  a bowthruster  is  to  be  integrated  into  the  flotilla,  its  location  and 
maximum  (zero-speed)  thrust  are  given  in  the  COEF  4 line  of  data.  The  thrust 
versus  speed  characteristics  are  given  in  two  COEF  5 lines,  in  terms  of  points 
describing  a curve  such  as  shown  in  Figure  1.  Note  that  the  ordinates  of  the 
curve  are  nomalized  by  the  zero  speed  thrust,  thus  different  bowthruster  rat- 
ings can  be  evaluated  by  changing  only  a single  line  of  data. 

The  integration  control  parameters,  given  in  the  INIT  1 data  set  specify 
the  time  intervals  of  the  interpolation.  The  step  size  determines  how 
frequently  the  solution  sumnary  is  printed  out  and  the  rudder  commands  are 
requested.  It  is  most  convenient  to  keep  the  step  size  as  a power  of  two, 
since  it  may  be  halved  successively  by  the  program  to  meet  the  error  con- 
trol limits.  The  values  used  in  the  Sample  Input  file,  Appendix  A,  are 
good  nominal  values  for  most  usage,  although  others  can  be  substituted  if 
desired. 

The  INIT  0 and  STER  data  sets  specify  the  initial  conditions  of  the  flo- 
tilla. All  values  are  self  evident  except  for  DRADIN  and  BOWCLR.  These 
two  parameters  define  the  initial  radial  offset,  dR,  of  the  barge  from  the 
center  line  of  the  channel  by  the  expression 

dR  «=  DRADIN  * (BOWCLR  * Ws  - (Ws  - WT)  / 2 ) 


BMAX 


Figure  1. 
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where 

Wg  = the  segment  width  at  its  beginning 

WT  = the  flotilla  width 

The  values  given  in  the  sample  input  data  file  will  serve  to  place  the 
flotilla  approximately  in  the  middle  of  the  channel  initially. 

Once  the  input  data  file  is  set  up,  the  simulation  program  can  be  ex- 
ecuted. The  program  will  read  all  data  for  the  first  case,  down  to  the 
TRIP  delimiter  line.  Execution  of  the  case  then  begins,  and  the  first  out- 
put is  simply  a printout  of  the  input  data.  This  can  be  directed  to  the 
users  terminal,  to  permit  checking  of  the  input  data,  or  if  the  user  is 
confident  that  his  input  file  is  correct,  the  listing  is  directed  to  an 
output  file,  to  be  listed  at  the  conclusion  of  the  run  with  a summary  of 
the  maneuver.  Samples  of  both  of  these  printouts  are  given  in  Appendix  B, 
along  with  a typical  display  that  is  generated  during  the  command  prompt. 

Rudder  and  propeller  commands  are  entered  in  response  to  the  prompts 
RUDDER  and  THROTTLE,  respectively.  The  steering  and  flanking  rudder  posi- 
tions (in  degrees)  are  entered  in  response  to  the  first  and  the  port,  star- 
board and  centerline  (as  applicable) .propeller  speed  (r.p.m.) .upper  limits 
are  entered  in  response  to  the  latter.  The  present  status  of  the  barge  mo- 
tion (its  r.p.m.,  rudder  position,  radial  velocity  in  feet  per  second  and 
yaw  rate  in  degree  per  second)  is  also  printed  out,  to  aid  in  selecting  an 
appropriate  command  setting. 

The  r.p.m.  is  printed  out  for  each  propeller,  port,  starboard  and  cen- 
terline respectively.  If  the  towboat  is  twin  screw,  the  centerline  propeller 
r.p.m.  will  be  shown  as  zero.  In  the  case  of  a single  screw  towboat,  the 
r.p.m.  will  be  displayed  by  the  label  'RPM-P'.  The  radial  velocity,  labelled 
'V-RAD',  represents  the  velocity  of  the  flotilla  outward  from  the  center  of 
the  curve.  The  relative  yaw  rate,  labelled  'V-YAW',  represents  the  rate  at 
which  the  flotilla  is  turning  in  relation  to  the  radially  directed  vector. 

The  rudder  positions,  'S-RUD'  and  'F-RUD’  indicate  the  present  position,  in 
degrees , or  the  steering  and  flanking  rudders,  respectively. 

Two  options  are  available  through  the  rudder  command.  If  the  steering 
rudder  position  is  greater  than  90°,  the  simulation  will  "back-track"  to 
the  preceedlng  time  step,  so  that  the  most  recent  command  can  be  altered. 

If  the  rudder  position  is  less  than  -90°,  the  simulation  will  "move"  the 
flotilla  back  to  the  beginning  of  the  present  segment,  and  reinitialize 
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Its  motion  to  the  values  that  were  applicable  at  that  time.  Thus,  if  a 
maneuver  is  represented  by  several  major  segments,  each  stage  can  be  re- 
petitively analyzed  before  proceeding  with  the  simulation,  without  re-run- 
ning  the  entire  maneuver  from  the  start.  This  is  a most  convenient  feature 
for  analyzing  a maneuver  that  requires,  say,  a sharp  turn  followed  by  a 

' I 

narrow  passage  and  another  bend.  Each  successvie  stage  can  be  exhaustively 
simulated  to  determine  the  best  strategy  for  approaching  the  following 
segment. 

The  simulation  thus  offers  a wide  range  of  flexibility  in  selecting 
different  barge  configurations  and  channel  descriptions.  The  input  data 
requirements,  both  from  the  data  file  and  from  the  interactive  command 
mode,  are  simple  and  convenient.  The  periodic  display  of  the  barge  posi- 
tion, and  interrogation  for  rudder  and  throttle  commands  enables  the  user 
to  develop  a "feel"  for  how  the  barge  responds  to  commands  and  environ- 
mental effects.  Thus,  the  simulation  offers  a beneficial  means  of  evaluat- 
ing causes  and  counter  measures  relevant  to  inland  waterway  bridge  collisions . 
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VICKSBURG 

PASSAGE  USING  THREE 

SEGMENTS. 

SAMPLE  CASE. 

BARG  1 

2 

70.0 

70.0 

BARG  2 

3 

574.5 

112.5 

283.0 

56.0 

140.0 

28.0 

BARG  3 

0 

10.4 

13.8 

0.88 

0.80 

BARG  4 

2 

500.0 

250.0 

10.0 

50.0 

50.0 

0.6 

0.0 

BARG  5 

-4 

9.5 

7.5 

0.85 

0.0 

0.0 

INIT  0 

2 

100.0 

9.0 

0.0 

0.0 

0.0 

O 

o 

1. 

IN  IT  1 

20 

8.0 

1.0 

0.01 

0.01 

STER  0 

0 

0.50 

0.06 

0.06 

0.001 

0.05 

4.0 

ROUT  0 

3 

ROUT  1 

1 

0.0 

5200. 

7800. 

ROUT  2 

1 

20.0 

5000. 

6800. 

ROUT  1 

2 

Q.O 

5000. 

6800. 

ROUT  2 

2 

15.0 

5900. 

6700  . 

ROUT  3 

2 

10.0 

6200. 

7000. 

ROUT  4 

2 

20.0 

6000. 

8200. 

ROUT  5 

2 

20.0 

8200. 

9400. 

ROUT  6 

2 

15.0 

8200. 

9300. 

ROUT  7 

2 

15.0 

8000. 

9300. 

ROUT  8 

2 

15.0 

7400. 

9000. 

ROUT  1 

3 

0.0 

7400. 

9000. 

ROUT  2 

3 

30.0 

6000. 

9000. 

COEF  1 

1 

-0.00058 

-0.0011 

0.0 

0.0 

COEF  2 

1 

-0.0033 

-0.0011 

0.0 

0.0 

COEF  l 

2 

0.0 

-0.0072 

0.00026 

-0.0070 

COEF  2 

2 

0.0 

-0.012 

-0.0081 

-0.016 

TRIP  2 
QUIT 


RUDDER= 

5.  5. 
THROTTLE= 
150.  150. 


3.2  PERCENT  THRU  SEGMENT  2 

I 

I 

I 

I 

I 

I 

I 


SPEED  6.1 
RPM-P130 .0 
RPM-S130 . 0 
-CL  0.0 
V-RAD  0.0 
V-YAW  0 . 1 
S-RUD  5.0 
F-RUD  5.0 


RUDDER= 

5.  5. 
THROTTLE= 
150.  150. 


B 

B 

B 

B 

B 

B 

T 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

0 

O 

O 

O 

0 

O 

O 

O 

O 

O 
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VICKSBURG  PASSAGE  USING  THREE  SEGMENTS.  SAMPLE  CASE. 

■ 


CHARACTERISTICS  OF  TOWBOAT  AND  3 LONG  BY  2 WIDE  BARGE  TOW 


TOW  . 

BOAT 

LENGTH  OVERALL 

575. 

113 

WIDTH 

70. 

70 

DRAFT 

10.40 

13 

BLOCK  COEFFICIENT 

0.880 

0 

L C G (FORWARD) 

283.0 

56 

GYRADIUS 

140.0 

28 

PROPULSION  AND  RUDDER  CHARACTERISTICS  AND  STEERING  CRITERIA 


NUMBER  OF  SHAFTS  2 

HORSEPOWER  PER  SHAFT  500. 

MAXIMUM  RPM  250. 

SHAFT  OFFSET  FROM  CL  10. 

BLADES  PER  PROPELLER  -4 

DIAMETER  9.5 

PITCH  7.5 

AREA  RATIO  0.850 

WAKE  FRACTION  0.0000 

THRUST  DEDUCTION  0.0000 

AREA  PER  STEERING  RUDDER  50. 

AREA  PER  FLANKING  RUDDER  50. 

MAXIMUM  RUDDER  ANGLE  0.600 

MAXIMUM  RUDDER  RATE  0.060 


INITIAL  BARGE/TOWBOAT  VELOCITY  AND  ORIENTATION 


BOWCLR  SPEED  GAMMA  YAW  YAW  RADIAL 

ANGLE  RATE  OFFSET 
0.50  9.00  0.0000  0.0000  0.0000  1.000 


INTEGRATION  CONTROL  PARAMETERS 


NCUTS  FIRST-STEP  STEP-SIZE  REL-ERROR  ABS-ERROR 

20  1.0000  8.0000  0.01000000  0.01000000 


DkOvOi/^tf^tnu^tf^iPi/'iippnpAPAPArA— 

0 4»^'fl'C^O'fi'fi'fl'C®®0&flD0)a'a'0'Cr(7‘^vC»C<C»C»C»C\C\0^v0vfi>flvtf\fivC>0»e>0\Ov0'fi>fiO^»fi\Cv0yfivC 
JCCCO(OOOOOCDtta)flDOOa3®0»®CD’r»^»V*CD®CDflDODttaO€DflDflDflDCOCOa)a)flD  OD'OO  ®flCflD®»®flDCDCDOOflD 


~nnnnnr»-hr*r-r^MnnnnintOif  iTu^ninnnN'-*- • •—*-*  — 0'a'C'C'^'r''HAPAHAPA* 
C (NiN(NtNfNi/'iTint/,'iri  pa  pa  pa  pa  pa  o O O O o pa  pa  pa  ^a  r-~  ua  y~i  ^iTiPr“f>,r'fv'r'fs'r'r'fx 


r-  O'  O'  O'  O'  ^ ^ ^ ^ ^ 


CJ 


-ooooococ3CDOcDininintninNiNNf'MNnH«riHnpnnnrt«««rtr«r»«»M«»,'«««««n« 

r^r»r,'r-r'r^r,->oovovONor*r*p-r,«p-ooaccooDcoO'0'0'0\C'0'0'0'0'0'0*o\0'0'0'0*0'0'0‘0'0'0'0'C>0' 

XiTurMrmin>H«--r-cDcca300flD»<r'rv»o'0'0'0'a'nnnnn»vvv»fNNrjfNr*r«rwvvvrw 

I I | | | I I I I I I I I I I i I I I I I I > I I 


-<flNO'CO'a''w'*'r'^nif>irii/UrirCDCD®COCDNNNCMNfNfMrj(N<NOOOOO(N(NfMWNr'M 

Or-^r-(ro^'r«0'0'0'0'0*OOOCO^-^^^-^r*r^r*r*r>*O^CO^£)vOV>w^^r’»vOvO^O^O'0^•^, 

CJ(N(N(N(NtM—  — — F-—  FI  — — 


r*  r~  NO  *£>  v£5  VC  vO 

1 ^ ^ r>  r*>  r*»  r*  r*> 

i pa  pa  vO  sO  vO  ^ 


rooooooooooioinirinrf'irunif>(ririnininininintririni^irininintr>i^ifiirurifMntr>iflifMO»o»nt/'inir 

a ^ rA  rA  pa  pa  n n n pa  n HA  ha  ha  fA  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  pa  ha  pa  pa  pa  pa  ha 

<fArApAPAfAPAPArAPAPA4nmtf^»n»nmininir.  *noooooooooo»r>ir»ininir>pAHAnPAHnpAPAPAPApAu^u"»ininir» 

g_r»r.r,r-.r^n*r*r'r*r^v^^^^^^,»^^©ooooooooo'3,^v^r,»f''r,^p*r»r'-r>*r'*r*r^r^^vv^v 
JoocoD2ooa)ceDor'r'r^f,*fsr^r^r,*r*r*oooooooooof'r^fs^‘f^ocDCDOODCDoocDoor^^‘r^r'r» 
UOOGOOCIOOOO*— ■»-*•-«»■■«.— •~H.—».—«~--0000000000— »-»oooooooooo*— —**—  —-r— 


oooooooooooooooooooocoooooooo< 

I I I I I I I I I I 


(OOOOOOOOOOOOOOOOOOOO 


r-r'0'ca'(Niror(NMF-rHr'rHv^M^p*iNf'vointN»ou^HnM(MPAi/,ivOFHnrH(7'r'r'aNn^p'ir'fNAA^ 

<(NA'«p(NccAipir'^F-icDin»ir>or'f,-*Oir>PAtN(NPAv»nkninir>«rPAovoo(NfMifivovovf'jff'ir'A>i)or^^O'r'(r 

XfNvi>0CDC''-t(NA^rinfAfN>-iOOcr'0NOOAninr'0'APininr'0'A^v00'Hn<rLnvflrvCDCD0'OOAMfN(\fAn 

<oooooAFHMH1Hl-HMHHOOHFiHHAAHiNfNNfSNnpopnn<r^Tw»r^^ininifurininiriir 

O 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

r~F-»r^-F-,cj''or'4f^F-<mr-«minf-<ocoir>vovo^,rMvo— <r^-^-<ir>a''^r^Aa3fM^HaDvDao'vir>04r«ivor^^ra'oa'r>4r'CAjpAaD 
<rMO'v£>^roo*-A^*-p*®in«»,»-tcoa,>r«iH®pvmooa'so^^Hoovo^Hco^N£)^AH(Nor>-HAoopAeD<Nr*ooo^-»-«o 

(HH(Nv^>coo'«-'fNPAvAHa'vx>^(NOAninvflvOi/iiPirmvvvc(NHOfNvinva\or^r^cDO^a'OAfMn^in 
UJOOOOOO*— i^HOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO— — 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I I I I I I I I I f I I I I I I I I I I I I I I I I I I I I I I I I I 

• 

ZDv£>aop^cDuAPAPAi(A|PAmvor^r^vo,o-fMOCDf^vor,,^0'PAO>vovop^oir>*— ‘OOHAoo^rc«<r^— 

MP'invviT'0r'CD(riOHfNrAcinu>Ar'ffia'O»-<nvyccDOHiincDr-vr'O^ffl(Nv0O»cDfNr'«vDOvovcD 

Dj'O'CMy'ff'CJ'a'ff'O'OCOOOOOOOOOHrt-.HHHININMfNPArAnv^rifunvO^VOAMEaDCMT'a^OO 

<LOvr>cninuAir>ir>uAuAu3vou?vovc>»ovc>vDvo\£>vDvovovovo\i3vovovosov^vDvDvovovcvovovovov^N£>vovovovo^v£>vor^r~ 

oc 

nr'voo(NrHr^(Nr'aMTMNr'VO'0'fN^r(j\0'a3a>TnAinofNHvconPAinaD^,v£fNO'CO'»roa3PAPAn\o(No 

QH(Ar*cDNO'co<T'fNaoo'ff>voo>o>CDvOHPANONnpApnAinvo(M^HHHOiriricDOov(Mr«An»0'0'®^co 
<oca'r'«ccoA(NNO'r^in^r^(NOHirtHOo^yco>f«j^»ccDO'00'f^',,,i\orHn^TfNor'n®fMor,*p,nocD 
ocfM^»r»r»a'OrMPA^^pA<Nr-«ooooo»-icN^vocoorAtr>r^a'H-«^\oo'<N^^«<7NO*— •iNPAPn«r^»n\ovor^r^r^ 


>OCOOO*HAAHHAHHAAAAAHAr^A( 


• NCNNNnnpAn^'w^»inif>»n»ninin»oir'»n»n>n»njr< 


pA\o^’pnovOHArMHAr»o%ovo\©^vo<T»oo-HNO*np^>r*PAPAo^rvo^,o^ir>r^flD'*,>p*®0'p»HOu^m«--<NivO'0*-«^o 
*-<c7''^r-ccacvo^-^r^rMPAOHNr^in^rmcDPAQOfM»A<>4r-ta'r'^^comaor^HAr^oH^ir>oo^-^sor~eDa'C'r'PAoo-^pA 
(/)oHHHriHHHOooAfsninr'ff'A<ftf(T>HA^>ocDON'^irui>MDcoa'(J'0'ff'Oooooooooo'a'tt 
Q*00000000000000000>— ••A^r-*rMf\J(NfNfMPAPAPAPArAPAPAPAPAPAPA/A^^^«r^5C-^^r^P»'PArA 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOl 


>0000000000000 


PnnAvnPAPAnnfAPAfMNN<N(N(NOJ(M(NN(NrMN(M(MMNfMNfMPOPA(AnPnnAPnPnnnvVirPAPnpAPAPS 

—>OOOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOOOOOOOO 

C I I I I I I I I t I I I I I I I I I I I I I I I I I I I I I I I I I I I t I I I I I 1 I 1 1 I I I I 

<U)UUUUU.'UUb:UUUb)UUUUUUUUUUUUb}UUUUUUUUUUUUUUU(L'UUUUUUUU 

WPA^o>vcA(Nff'O'OiAr*»cH»('jifl<oin(SMnoMOO'nff>>fiNO'»O'OP*'r,''OHC0r'(N«(T'A»(Ntrnr'vtrinA 

nr*»f-*F»rA  — pAfAPApA<Mr^<Nr^oHAmr»a'FA<Nr^^HA<MrM— iHOO'AiinnoH»oa\hNvOHino>\cff'hO'Off>o 

^flD«r(NPArt(MPnvir«OinHHM(MNNNfAPAr4NN(NNN(M(M(NNHO'»flVNnPnnPAOMNAO'»AP>HOM»A 



OOOOOOGOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I I I I I I I I I I I I 

<rNknAOHurxAvr>ovr>0'»r'Air(y\(NvoO'rflDfNinavpAvoo^voc^(Ninr'OfN»'coDO<Ni*'invooDC'o  — 
~ — rMrAinsor-aoo— KNviovocoooN'Av^MJ'OAM^inhcooA i<NAin»or'<y'0- <rNPA»nvcr~coiT>o-—  a v 
aoooooooAHHHHHHHNN(NMNNNn(AfAPA(APAA»*vr’rwtiftinininininiOrfMn»fi>ovo»c 



•<00000000000000000000000000000000000000000000000000 

MJNNfN(NfMNfXWrMr'irJ(NM<NNM(N(N(N(MrJM(NMr«i<NPN(N(Nr4(Nir4fNNN(N(N{N(NfNNfNfMPNr>HNfNPNM 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
^ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

Oh.'b.'UUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUU(iHL'UU 
Z^i/^r^-aoo— p>JfMrA^rinvor^ooa'0^r^PAPA^rPAPAPArN^-«F-^oor>inpA— •0'r>vr>^rHAr^oa'ODvr>uA^rPA^oc7'ao 
<in»nininvo\ovcsoifi^»o»ovo»o*flr'f'r'f^r*r*r'r^r'r'P'r»r^\o\o«o»ovo>fl»runinir>»/'»A'^‘»»»,»^,r»'An 



00000000000000000000000000000000000000000000000000 

li.'eD'0,»NC®vO«NO<D«)VNO«'OfNOeDvCTNOaMfl*NO(D>e«fNOflDvC*NOCD\C'r(NCflO»OrNO 
H A(MAwi/>vo^®eoo>o^(N(N  r^^vovcr>®(roOANPA»viOyfir>a5®<y'OAf'4fNAv^v£v£ir^a30'c 
A-  ^^HpNHH»-MHF-HHfMNNNNNNNNNNN(NPnfAnP^PnnnnAAAr>r 

t~ 

D,F-<(Nn»»r»\cr'®(5'OAPNPn»if\®p'®(T'OA(NAV1nxer“®oOF-(NAv»r®r'®C'0--HgpA^^\or>a3C'0 
U HFH^HAF-p*AAF-r4NNNNNiNNNNPnpAnpAnnpA(AfAn«r'rw«w»wir 

f- 


STEP  TIME  VIANG)  ALPHA  V (HEAD)  PSI  V ( RAD)  RADIUS  BETA  GAMA  DELTA  RPM  CLIN  XLIN  CLOUT  XLOUT 

SI  408 . 0. 135E-02  0.6S2  -0.508E-03  0.3780  5.6244  7133.  -0.1517  0.5297  0.0000  80.  977.  -491.  378.  1865. 
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BARMAN*** 

BARGE  MANEUVERING  SIMULATION 
***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN  , LPRT , MOJ4ITR 

COMMON/CMCNTR/NCARD ,NCASE ,NERORS ,NPAGE , KPATH ,KPRINT , LINECT 

CALL  BEGIN 

CONTINUE 

CALL  INPUT 

CALL  INITAL 

CALL  ROUTCK 

CALL  TRPINT 

CALL  TRPOUT 

GO  TO  10 

END 

SUBROUTINE  BEGIN 

*** IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT , MONITR 

COMMON/CMCNTR/NCARD , NCASE , NERORS , NPAGE , KPATH , KPRINT , LINECT 
COMMON/CMDATA/TEMP ( 7) ,NHEAD ( 20 ) , LABEL, ITAG, JTAG 
COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD , NPROP , NBLAD 
COMMON/CMINIT/SPDIN ,CURENT , GAMIN , HEADIN , CDOTIN , DRADIN , INSEG 
COMMON/CMROUT/ SDAT (11,10,2) ,CUR ( 11 , 10 , 3) ,SANG(11,10) ,NANG(10) , NSEG 
COMMON/CMCNST/RHO ,GRAV , PI 

COMMON/CMCOEF/A ( 10 ) ,B (10)  , AA (10)  ,BB(10)  ,COEF(5) 

COMMON /CMC LOG/NPOLY , NSTRI P , KORPOL , KORSTR 

COMMON/ CMBOWT/BTHRUS ( 7 ) ,BTSPD ( 7 ) ,BTMAX , BTPOS , BTGAIN , NBTSPD , IDBT 
LOGICAL  NPOLY , NSTRI P , KORPOL , KORSTR , IDBT 
THIS  SUBROUTINE  ASSIGNS  INITIAL  VALUES  TO  SEVERAL  CONSTANTS 
IN=20 

WRITE (5 ,88) 

FORMAT ('  ENTER  INPUT  FILE  NAME') 

READ (5,89) BARIN 
FORMAT (A5) 

CALL  IFILE (20, BARIN) 

CALL  OFILE(21, 'BAROUT') 

LPRT=21 

MONITR=5 

RHO=l . 99 

GRAV=32 . 2 

PI=3. 141592 

READ  (IN, 1000)  NHEAD 

KORPOL=. FALSE. 

KORSTR=. FALSE. 

IDBT=. FALSE. 

NCARD=1 

NCASE=0 

NPAGE=0 

BNUL=0 . 0 

WFRAC=BNUL 

TDDUC=BNUL 
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10500  CURENT=BNUL 

10510  GAMIN=BNUL 

10520  HEADIN=BNUL 

10530  CDOTIN=BNUL 

10540  DRADIN=BNUL 

10550  DO  10  J=l,10 

10560  NANG ( J ) =0 

10570  DO  10  1=1,11 

10580  SANG (I , J) =BNUL 

10590  DO  10  K=1 , 2 

10600  SDAT ( I , J , K ) =BNUL 

10610  10  CONTINUE 

10620  RETURN 

10630  1000  FORMAT  (20A4) 

10640  END 

10650  SUBROUTINE  INPUT 

10660C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

10670  COMMON/CMPRNT/IN , LPRT , MONITR 

10680  COMMON/CMCNTR/NCARD , NCASE , NERORS , NPAGE , KPATH , KPRINT , LINECT 

10690  COMMON/CMDATA/TEMP ( 7 ) ,NHEAD (20) , LAB EL, IT AG, JT AG 

10700  COMMON/ CMWIDE/TOWWID , BOTWID , NBARW 

10710  COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 

10720  COMMON/CMCHAR/TOWDRF , BOTDRF , TOWBC , BOTBC 

10730  COMMON/CM PROP/SHP , DPROP , PI TCH , ARAT , WFRAC , TDDUC , PD , NPROP , NBLAD 

10740  COMMON/CMINIT/SPDIN, CURENT, GAMIN, HEADIN, CDOTIN, DRADIN, INSEG 

10750  COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 

10760  COMMON/ CMWIND/ SWIND (10) ,DWIND(10) 

10770  COMMON/CMSTAT/RPM ( 3 ) , SPEED , DELTA ( 2 , 3 ) , GAMMA , HEAD ,CDOT, DRAD , 

10780  &RPMHI (3) ,ISEG 

10790  COMMON/CMNTIG/FIRSTP , STEP , EPS , AB , NCUTS 

10800  COMMON/CMRUDR/ARUD (2,3) , DELMAX , DLDTMX , OFSET 

10810  COMMON/CMSTER/BOWCLR, STRSLO , STRSLI , STRCOR , STRBK , RGAIN , RPMAX 

10820  COMMON/CMCNST/RHO ,GRAV , PI 

10830  COMMON/CMCOEF/A ( 10 ) ,B(10) ,AA(10) ,BB(10) ,COEF(5) 

10840  COMMON/ CMCLOG/NPOLY , NSTRI P , KORPOL , KORSTR 

10850  COMMON/CMBOWT/BTHRUS ( 7 ) ,BTSPD(7) ,BTMAX , BTPOS , BTGAIN , NBTSPD , IDBT 

10860  LOGICAL  NPOLY ,NSTRIP , KORPOL, KORSTR, IDBT 

10870  DIMENSION  NAME ( 9 ) 

10880  DATA  NAME/4HBARG , 4HINIT , 4HROUT , 4HNTRP , 4HCOEF , 4HPARA , 4HSTER , 

10890  &4HTRIP , 4HQUIT/ 

10900C***  THIS  SUBROUTINE  READS  A COMPLETE  SET  OF  INPUT  DATA  ON  EACH  PASS 
10910  1 CONTINUE 

10920  NCASE=NCASE+1 

10930  NPOLY=. FALSE. 

10940  NSTRIP=. FALSE. 

10950  NERORS=0 

10960  10  READ  (IN, 1000)  LABEL, ITAG , JTAG , (TEMP ( I ) , 1=1 , 7 ) 

10970  NCARD=NCARD+1 

10980  JUMP=10 

10990  DO  20  1=1,9 
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11000  20  IF  (LABEL  .EQ.  NAME ( I ) ) JUMP=I 

11010  GOTO  (100,200,300,400,500,600,700,800,900,30) , JUMP 

11020  30  CALL  ERSTOP  (1) 

11030  100  CONTINUE 

11040C***  PRESENT  DATA  CARD  IS  A BARGE/TOW  DATA  CARD 
11050  IF  (ITAG  .LT.  1 .OR.  ITAG  .GT.  5)  CALL  ERSTOP  (2) 

11060  GOTO  (110,120,130,140,150)  , ITAG 

11070  110  NBARW=JTAG 

11080  TOWWID=TEMP ( 1 ) 

11090  BOTWID=TEMP ( 2 ) 

11100  GO  TO  10 

11110  120  NBARL=JTAG 

11120  TOWLEN=TEMP ( 1 ) 

11130  BOTLEN=TEMP ( 2 ) 

11140  CGTOW=TEMP ( 3 ) 

11150  CGBOT=TEMP ( 4 ) 

11160  TOWK=TEMP (5) 

11170  BOTK=TEMP (6) 

11180  GO  TO  10 

11190  130  TOWDRF=TEMP ( 1 ) 

11200  BOTDRF=TEMP ( 2 ) 

11210  TOWBC=TEMP ( 3 ) 

11220  BOTBC=TEMP ( 4 ) 

11230  GO  TO  10 

11240  140  NPROP=JTAG 

11250  SHP=TEMP ( 1 ) 

11260  RPMAX=TEMP ( 2 ) 

11270  OFSET=TEMP ( 3 ) 

11280  DO  145  1=1,3 

11290  ARUD(1 , I) =TEMP(4) 

11300  145  ARUD (2,1) =TEMP ( 5 ) 

11310  DELMAX=TEMP ( 6 ) 

11320  DLDTMX=TEMP (7 ) 

11330  GO  TO  10 

11340  150  NBLAD=JTAG 

11350  DPROP=TEMP ( 1 ) 

11360  PITCH=TEMP (2) 

11370  ARAT=TEMP (3) 

11380  WFR AC=TEMP ( 4 ) 

11390  TDDUC=TEMP ( 5 ) 

11400  GO  TO  10 

11410  200  CONTINUE 

11420C***  PRESENT  DATA  CARD  GIVES  INITIAL  CONDITIONS  FOR  THIS  CASE 
11430  IF  (ITAG  .GT.  0)  GO  TO  210 

11440  INSEG=JTAG 

11450  SPDIN=TEMP (2) 

11460  GAMIN=TEMP ( 4 ) 

11470  HEADIN=TEMP(5) 

11480  CDOTIN=TEMP (6) 

11490  DRADIN=TEMP ( 7 ) 
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11500 

GO  TO  10 

11510 

210 

NCUTS=JTAG 

11520 

STEP=TEMP(1) 

11530 

FIRSTP=TEMP (2) 

11540 

EPS=TEMP ( 3 ) 

11550 

AB=TEMP ( 4 ) 

11560 

GO  TO  10 

11570 

300 

CONTINUE 

11580C 

*** 

PRESENT  DATA  CARD  DESCRIBES  CHANNEL  CHARACTERISTICS 

11590 

IF  (ITAG  .GT.  0)  GO  TO  310 

11600 

NSEG=JTAG 

11610 

GO  TO  10 

11620 

310 

CONTINUE 

11630 

SANG ( ITAG , JTAG) =TEMP ( 1 ) 

11640 

SDAT ( ITAG , JTAG , 1 ) =TEMP ( 2 ) 

11650 

SDAT ( ITAG , JTAG , 2 ) =TEMP ( 3 ) 

11660 

CUR ( ITAG , JTAG , 1 ) =TEMP ( 4 ) 

11670 

CUR (ITAG, JTAG, 2)  =TEMP(5) 

11680 

CUR ( ITAG , JTAG , 3 ) =TEMP ( 6 ) 

11690 

I F ( ITAG . NE . 1 ) GOTO  10 

11700 

SWIND ( JTAG) =TEMP(1) *1.47 

11710 

DWIND (JTAG) =TEMP ( 7 ) 

11720 

400 

CONTINUE 

11730 

GO  TO  10 

11740 

500 

CONTINUE 

11750C 

*** 

PRESENT  DATA  CARD  HAS  HYDRODYNAMIC 

COEFFICIENTS 

11760 

IF  (ITAG  .LT.  1 .OR.  ITAG  .GT.  5) 

CALL  ERSTOP  (2) 

11770 

GO  TO  (510,520,530,540,550),  ITAG 

11780 

510 

CONTINUE 

11790 

KORPOL= . TRUE . 

11800 

NPOLY® .TRUE. 

• 

11810 

IF  (NSTRIP)  CALL  ERSTOP  (3) 

11820 

IF  (JTAG  .GT.  1)  GO  TO  515 

11830 

AA (2) =TEMP ( 1 ) 

11840 

AA ( 3 ) =TEMP (2) 

11850 

AA (9) =TEMP (3) 

11860 

AA (10) *TEMP ( 4 ) 

11870 

GO  TO  10 

11880 

515 

CONTINUE 

11890 

AA ( 1 ) =TEMP ( 1 ) 

11900 

AA ( 5) =TEMP (2) 

11910 

' 

AA (6 ) =TEMP (3 ) 

11920 

AA ( 7 ) =TEMP (4) 

11930 

GO  TO  10 

11940 

520 

CONTINUE 

11950 

KORPOL® . TRUE . 

11960 

NPOLY=.TRUE. 

11970 

IF  (NSTRIP)  CALL  ERSTOP  (3) 

11980 

IF  (JTAG  .GT.  1)  GO  TO  525 

11990 

BB (2) =TEMP(1) 
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12000  BB (3) =TEMP (2) 

12010  BB ( 9 ) =TEMP ( 3 ) 

12020  BB ( 10 ) =TEMP ( 4 ) 

12030  GO  TO  10 

12040  525  CONTINUE 

12050  BB ( 1 ) =TEMP ( 1 ) 

12060  BB ( 5 ) =TEMP ( 2 ) 

12070  ' BB ( 6 ) =TEMP ( 3 ) 

12080  BB (7) =TEMP (4) 

12090  GO  TO  10 

12100  530  CONTINUE 

12110  KORSTR= . TRUE . 

12120  N STRIP* .TRUE . 

12130  IF  (NPOLY)  CALL  ERSTOP  (3) 

12140  DO  535  1=1,5 

12150  535  COEF ( I ) =TEMP ( I ) 

12160  IF  ( JTAG  .NE.  10)  GO  TO  10 

12170  NSTRIP=. FALSE. 

12180  KORPOL=. FALSE. 

12190  KORSTR=. FALSE. 

12200  GO  TO  10 

12210  540  CONTINUE 

12220  IF  (JTAG  .GT.  O)  GO  TO  545 

12230  IDBT=. FALSE. 

12240  GO  TO  10 

12250  545  CONTINUE 

12260  IDBT=.TRUE. 

12270  NBTSPD=JTAG 

12280  BTMAX=TEMP(1) 

12290  BTPOS=TEMP ( 2 ) 

12300  BTGAIN=TEMP (3) 

12310  GO  TO  10 

12320  550  CONTINUE 

12330  IF  (JTAG  .GT.  1)  GO  TO  555 

12340  DO  551  1=1,7 

12350  551  BTSPD ( I ) =TEMP ( I ) 

12360  GO  TO  10 

12370  555  CONTINUE 

12380  DO  556  1=1,7 

12390  556  BTHRUS ( I ) =TEMP ( I ) 

12400  GO  TO  10 

12410  600  CONTINUE 

12420  GO  TO  10 

12430  700  CONTINUE 

12440  BOWCLR=TEMP ( 1 ) 

12450  GO  TO  10 

12460  800  CONTINUE 

12470  K PATH* I TAG 

12480  KPRINT=JTAG 

12490  RETURN 


BARGE  PAGE  6 


12500  900  STOP 

12510  1000  FORMAT  ( A4 , IX , 11 , 1 3 , IX , 7F10 . 5 ) 

12520  END 

12530  SUBROUTINE  INITAL 

12540C***  ***IMPLICIT  REAL* 8 (A-H,0-Z) 

12550  COMMON/CMPRNT/IN , LPRT , MON ITR 

12560  COMMON/ CMCNTR/NCARD , NCASE , NERORS , NPAGE , KPATH , KPRINT , LINECT 

12570  COMMON /CM DAT A/TEMP (7) ,NHEAD ( 20 ) , LABEL, IT AG , JTAG 

12580  COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

12590  COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 

12600  COMMON/CMCHAR/TOWDRF , BOTDRF , TOWBC , BOTBC 

12610  COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD , NPROP , NBLAD 

12620  COMMON/CMINIT/SPDIN ,CURENT, GAMIN ,HEADIN ,CDOTIN ,DRADIN , INSEG 

12630  COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) , NSEG 

12640  COMMON/CMCNST/RHO ,GRAV , PI 

12650  COMMON/CMDISP/TOWDSP,BOTDSP ,GYRAD ,TMASS , ZNERTA , TLEN , TLEN2 , TLEN3 

12660  COMMON/ CMSTAT/RPM ( 3 ) , SPEED , DELTA ( 2 , 3 ) , GAMMA, HEAD, CDOT,DRAD, 

12670  &RPMHI  (3 ) , ISEG 

12680  COMMON/CMSTER/BOWCLR, STRSLO, STRSLI , STRCOR , STRBK , RGAIN , RPMAX 

12690  COMMON/CMSIZ E/CGAFT , CGFWD , EFLEN , EFDRF , EFBEAM , TONS 

12700C***  THIS  SUBROUTINE  ASSIGNS  INITIAL  VALUES  TO  SEVERAL  VARIABLES 
12710  TOWDSP=TOWLEN*TOWWID*TOWDRF*TOWBC 

12720  BOTDSP=BOTLEN*BOTWID*BOTDRF* BOTBC 

12730  TMASS= (TOWDSP+BOTDSP) *RHO 

12740  CGAFT= (TOWDSP*CGTOW-BOTDSP* (BOTLEN-CGBOT) ) / (TOWDSP+BOTDSP) 

12750  CGFWD=TOWLEN-CGAFT 

12760  PD=PITCH/DPROP 

12770  ZNERTA=BOTDSP* (BOTK*BOTK+ (CGAFT+BOTLEN-CGBOT) **2) 

12780  ZNERTA=ZNERTA+TOWDSP* (TOWK*TOWK+ (CGTOW-CGAFT) **2) 

12790  ZNERTA= ZNERTA* RHO 

12800  EFLEN=TOWLEN+BOTLEN*BOTWID/TOWWID 

12810  EFDRF=TOWDRF 

12820  EFBEAM=TOWWID 

12830  TONS= (TOWDSP+BOTDSP) *RHO*GRAV/2000 . 0 

12840  CALL  DATOUT 

12850  CALL  HYCOEF 

12860  DO  5 J=1 ,NSEG 

12870  SANG (1 , J) =0 . 

12880  NANG ( J) =0 

12890  DO  5 1=2,11 

12900  I F ( SANG ( I , J ) . EQ .0.0) GOTO 5 

12910  SANG ( 1 , J ) =SANG ( 1 , J ) +S ANG ( I , J ) 

12920  NANG ( J) =NANG ( J) +1 

12930  5 CONTINUE 

12940  SPEED=SPDIN 

12950  GAMMA=GAMIN 

12960  HEAD=HEADIN 

12970  CDOT=CDOTIN 

12980  I-INSEG-1 

12990  10  CONTINUE 
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13000  1=1+1 

13010  IF  (I  .GT.  NS EG)  GO  TO  100 

13020  SEGW=SDAT (1,1,2) -SDAT (1,1,1) 

13030  SEGD= (SDAT (1,1,1) +SDAT ( 1 , I , 2) ) /2 . 

13040  DRAD=DRADIN  * (BOWCLR+SEGW- (SEGW-TOWWID) /2.0) 

13050  RADl=SEGD+DRAD 

13060  ALDOT=SPEED/RADl 

13070  IF  (SANG (1 , INSEG)  .GT.  0.0)  GO  TO  100 

13080  ALDOT=- ALDOT 

13090  GAMMA=GAMMA+PI 

13100  HEAD=HEAD+PI 

13110  100  CONTINUE 

13120  CDOT=CDOT+ALDOT 

13130  RETURN 

13140  END 

13150  SUBROUTINE  ROUTCK 

13160C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

13170  COMMON/CM PRNT/IN , LPRT , MONITR 

13180  COMMON/CMINIT/SPDIN , CURENT , GAMIN , HEADIN , CDOTIN , DRADIN , INSEG 

13190  COMMON/ CMROUT/ SDAT (11,10,2) ,CUR ( 11 , 10 , 3 ) ,SANG(11,10) ,NANG(10) , NSEG 

13200  COMMON/CMCNST/RHO ,GRAV , PI 

13210  COMMON /CMTRIP/TRPDST , ENDHED , STRDST , RSLTNT 

13220C***  THIS  SUBROUTINE  COMPUTES  THE  NET  DISTANCE  AND  HEADING  CHANGE 

13230C***  DURING  A COMPLETE  TRIP 

13240  XBASE=0 . 0 

13250  YBASE=0 . 0 

13260  TRPANG=0 . 0 

13270  TRPDST=0 . 0 

13280  DO  100  J=INSEG ,NSEG 

13290  NSUB=NANG ( J) +1 

13300  DO  100  1=2 ,NSUB 

13310  ANGRAD=SANG(I ,J)/57.3 

13320  SEGD= ( SDAT ( I , J , 1 ) +SDAT ( I , J , 2 ) ) /2 . 

13330  TRPDST=TRPDST+SEGD*ABS (ANGRAD) 

13340  SGNR=SEGD*ANGRAD/ABS (ANGRAD) 

13350  XPRIM=SGNR*SIN (ANGRAD) 

13360  YPRIM=SGNR* (1.0-COS (ANGRAD) ) 

•13370  XBASE=XPRIM*COS (TRPANG) -YPRIM+SIN (TRPANG) +XBASE 

13380  YBASE=YPRIM*COS (TRPANG) +XPRIM*SIN (TRPANG) +YBASE 

13390  30  CONTINUE 

.13  400  TRPANG=TRPANG+ ANGRAD 

13410  100  CONTINUE 

13420  ENDHED=TRPANG*180.0/PI 

13430  STRDST=SQRT ( XBASE *XBASE+YBASE*YBASE) 

13440  RSLTNT=ATAN2 (YBASE , XBASE) *180. 0/PI 

13450C***  THIS  IS  A TEMPORARY  PRINT  STATEMENT 
13460  WRITE  (LPRT, 1000)  TRPDST, ENDHED , STRDST , RSLTNT 

13470  1000  FORMAT  ('  CURVE , HEAD, STRAIGHT, NET= ', 4F15 . 5) 

13480  RETURN 

13490  END 
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13500 

13510C*** 

13520 

13530 

13540 

13550 

13560 

13570 

13580 

13590 

13600 

13610C*** 

13620C*** 

13630 

13640 

13650 

13660 

13670  1000 
13680 
13690 
13700 
13710 
13720 
13730 
13740 
13750 
13760 
13770 
13780 
13790 
13800 
13810 
13820 
13830 
13840 
13850 
13860 
13870 
13880 
13890 
13900 
13910 
13920 
13930 
13940 
13950 
13960 
13970 
13980C*** 
13990 


SUBROUTINE  TRPINT 
***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT , MON ITR 

COMMON/CMCNTR/NC ARD , NC ASE , NERORS , NPAGE , KPATH , KPRINT , LINECT 
COMMON/CMINIT/SPDIN ,CURENT , GAMIN , HEADIN ,CDOTIN , DRADIN , INSEG 
COMMON/CMROUT/SDAT (11 ,10,2)  ,CUR (11 , 10 , 3)  ,SANG(11,10)  ,NANG(10)  , NSEG 
COMMON/CMSTAT/RPM (3) , SPEED , DELTA ( 2 , 3) , GAMMA , HEAD ,CDOT , DRAD , 

&RPMHI ( 3 ) , ISEG 

COMMON/CMTIME/ACUMT (100)  ,STIME (100)  , NXTCRV , LSTCRV,NXTREV, LSTREV 
COMMON/CMSCR1/ALRCOM 

LOGICAL  NXTCRV, LSTCRV,NXTREV, LSTREV 
THIS  SUBROUTINE  CONTROLS  THE  COMPUTATION  OF  ELAPSED  TIME 
AS  THE  PROGRAM  PROGRESSES  ALONG  A ROUTE 
TIMEL=0 . 0 
NSEGM=NSEG-1 
DO  100  ISEG=INSEG ,NSEGM 
WRITE (LPRT, 1000) ISEG 
FORMAT  (//'  INTEGRATION 


V(ANG) 

RADIUS 


XLIN  CLOUT  XLOUT ' ) 


GO  TO  20 


20 


30 


100 


&'  STEP  TIME 
&'  PSI  V(RAD) 

& ' RPM  CLIN 
CONTINUE 
ALRCOM=0 .0 
TSEG=0 . 0 
ALPHA=0 . 0 
NEMORY=0 
NXTCRV=. FALSE. 

LSTCRV=. FALSE. 

NXTREV=. FALSE. 

LSTREV=. FALSE. 

J=ISEG+1 
IF  (J  .GT.  NSEG) 

IF  (SANG (1,J)*SANG(1, ISEG) 

NXTCRV= . TRUE . 

CONTINUE 
J=ISEG-1 

IF  (J  .LT.  INSEG)  GO  TO  30 
IF  (SANG (1,J)*SANG(1, ISEG) 

LSTCRV= . TRUE . 

CONTINUE 

CALL  SEGINT  (TSEG , ALPHA, NEMORY) 

TIMEL=TIMEL+TSEG 

ACUMT (ISEG) =TIMEL 

STIME (ISEG) =TSEG 

CONTINUE 

RETURN 

END 

SUBROUTINE  SEGINT  (TSEG , ALPHA , NEMORY) 
***IMPLICIT  REAL* 8 (A-H,0-Z) 
COMMON/CMPRNT/IN , LPRT , MONITR 


THROUGH  SEGMENT  NUMBER', 1 5// 
ALPHA  V (HEAD)  ', 

BETA  GAMA 


DELTA  1 


LT.  0.0)  NXTREV= . TRUE . 


LT.  0.0)  LSTREV= . TRUE . 
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14000  COMMON/CMCNTR/NCARD ,NCASE ,NERORS ,NPAGE , KPATH , KPRINT , LINECT 

14010  COMMON/CM INIT/SPDIN ,CURENT , GAMIN , HEADIN , CDOTIN , DRADIN , INSEG 

14020  COMMON/CMROUT/ SDAT (11,10,2)  ,CUR ( 11 , 10 , 3 ) ,SANG(11,10)  ,NANG(10)  , NSEG 

14030  COMMON/CMCNST/RHO ,GRAV , PI 

14040  COMMON/CMPROP/DUM ( 7 ) ,NPROP,NBLAD 

14050  COMMON/CMDISP/TOWDSP ,BOTDSP ,GYRAD ,TMASS , ZNERTA ,TLEN , TLEN2 , TLEN3 

14060  COMMON/CMSTAT/RPM ( 3 ) , SPEED , DELTA ( 2 , 3) , GAMMA , HEAD ,CDOT , DRAD , 

14070  &RPMHI (3) , ISEG 

14080  COMMON/CMTRIP/TRPDST , ENDHED , STRDST , RSLTNT 

14090  COMMON/CMTIME/ACUMT (100)  ,STIME(100) ,NXTCRV , LSTCRV f NXTREV , LSTREV 

14100  COMMON/CMBEGN/BEGRPM (3) , BEGSPD , BEGDEL ( 2,3) ,BEGGAM , BEGHED ,BEGCDT , 

14110  SrBEGDRD 

14120  COMMON/CMSTER/BOWCLR, STRSLO , STRSLI , STRCOR , STRBK , RGAIN , RPMAX 

14130  COMMON/CMS I Z E/CG AFT , CGFWD , EFLEN , EFDRF , EFBE AM , TONS 

14140  COMMON/ CMCLER/ CLIN , XLIN , CLOUT , XLOUT , SIGN , CLBOW , BOWCL 

14150  COMMON/SCR2/NOUTS ( 8 ) ,POUTS(8) , NSTRS ,KNT 

14160  LOGICAL  NXTCRV, LSTCRV, NXTREV, LSTREV, FFALSE 

14170  DIMENSION  Y(6) 

14180C***  THIS  SUBROUTINE  CONTROLS  THE  INTEGRATION  OF  THE  EQUATIONS  OF 
14190C***  MOTION  THROUGH  ONE  SEGMENT  OF  A ROUTE 
14200  FFALSE=. FALSE. 

14210  DO  1 1=1 ,NPROP 

14220  RPMHI ( I ) =RPMAX 

14230  BEGDEL (1,1) =DELTA (1,1) 

14240  BEGDEL (2,1) =DELTA (2,1) 

14250  1 BEGRPM ( I ) =RPM ( I ) 

14260  BEGSPD=SPEED 

14270  BEGGAM=GAMMA 

14280  BEGHED=HEAD 

14290  BEGCDT=CDOT 

14300  BEGDRD=DRAD 

14310  10  CONTINUE 

14320  NSTRS=0 

14330  SIGN=SANG ( 1 , ISEG) /ABS (SANG ( 1 , ISEG) ) 

14340  CALL  SETY  (Y, ALPHA) 

14350  CALL  INTEG  (0 ,TSEG ,Y,J900 , FFALSE) 

14360  IF  (ALPHA  .EQ.  0.0)  NTIMES=0 

14370  20  CONTINUE 

14380  KNT=1 

14390  NTIMES=NTIMES+1 

14400  IF  (ABS (Y (2) ) .LT.  ABS (SANG ( 1 , ISEG) /57 . 3) ) GOTO  30 

14410  CALL  TIMCOR  (Y,TSEG, KPATH) 

14420  GO  TO  999 

14430  30  CONTINUE 

14440  IF ( ( 4+NTIMES) /5*5 .NE . 4+NTIMES) GOTO 31 

14450  29  CONTINUE 

14460  CALL  PLOTER ( Y (6 ) , Y (4 ) , Y ( 2 ) , ISEG) 

14470  118  CONTINUE 

14480  WRITE (5,311) 

14490  READ (22,*, ERR= 118) (DELTA (NR, 1) ,NR=1,2) 
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14500 

14510 

14520 

14530 

14540 

14550 

26 

14560 

119 

14570 

14580 

14590 

937 

14600 

14610 

14620 

31 

14630 

14640 

311 

14650 

312 

14660 

313 

14670 

14680 

14690 

32 

14700 

14710 

33 

14720 

14730 

14740 

14750 

14760 

14770 

14780 

14790 

14800 

14810 

14820 

35 

14830 

14840 

14850 

900 

14860 

14870 

910 

14880 

14890 

14900 

950 

14910 

999 

14920 

1000 

14930 

2000 

14940 

3000 

14950 

14960 

14970 

14980C 

*** 

14990 

WRITE (5,937) (DELTA (NR , 1 ) ,NR=1,2) 

DELTA (1,1) =DELTA (1,1) /57 . 3 
DELTA (2,1) =DELTA ( 2 , 1) /57 . 3 
DO  26  NR=2 ,NPROP 
DELTA ( 1 , NR) =DELTA ( 1 , 1 ) 

DELTA ( 2 ,NR) =DELTA (2,1) 

WRITE (5,313) 

READ ( 22 , * , ERR=119) (RPMHI (NR) ,NR=l,NPROP) 

WRITE (5,937) (RPMHI (NR) ,NR=l,NPROP) 

FORMAT (2F5.0) 

I F ( ABS ( DELTA (1,1) ) .GT. 1 . 57 ) G0T032 
CALL  SETIC  (Y ,TSEG, ALPHA, NTIMES ,NEMORY) 

CONTINUE 

CALL  INTEG  (6 ,TSEG , Y , J900 , FFALSE) 

FORMAT ( ’ RUDDER= ’ ) 

FORM AT (FI 0.5) 

FORMAT ( 1 THROTTLE = ' ) 

ALPHA=Y (2) 

GOTO 3 3 
CONTINUE 

CALL  SETIC  (Y,TSEG, ALPHA, NTIMES, NEMORY) 

CONTINUE 

GAMA=ATAN2 ( Y ( 5 ) ,Y(6)*Y(1)) 

IF  (SIGN  .LT.  0.0)  GAMA=GAMA+PI 

BETA=Y (4) -GAMA 

RPMAV= ( RPM ( 1 ) +RPM ( 2 ) ) /2 . 

DELTAS= (DELTA (1,1) +DELTA ( 1 , 2 ) ) /2 . 

WRITE  (LPRT, 2000)  NTIMES ,TSEG , ( Y ( I ) , 1=1 , 6) ,BETA , GAMA , 

&DELTAS , RPMAV , CLIN , XLIN , CLOUT , XLOUT 
IF (ABS (DELTA (1, 1) ) .GT . 1 . 57 ) G0T029 
IF  (NTIMES  .LT.  200)  GO  TO  20 
GOTO  (999,910,35) ,KPATH 
CONTINUE 
TSEG=-1 . 0 
GO  TO  950 
WRITE  (LPRT, 1000) 

GO  TO  999 
CONTINUE 
TSEG=0 .0 

CALL  TIMCOR  (Y,TSEG,4) 

WRITE  (LPRT, 3000) 

RETURN 

FORMAT  (’  ITERATION  FAILURE’) 

FORMAT  (I5,F7.0,E12.3,F8.3,E12.3,2F9.4,F8.0,3F9.4,5F7.0) 

FORMAT  ('  PRECEEDING  SEGMENT  NOT  COMPLETED  WITHIN  ALLOWED  LIMITS 
&) 

END 

SUBROUTINE  ERSTOP  (NGOOF) 

★★♦IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/ CMPRNT/IN , LPRT , MONITR 


I 
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15000  COMMON/CMCNTR/NCARD,NCASE ,NERORS ,NPAGE , KPATH , KPRINT , LINECT 

15010  COMMON/CMDATA/TEMP ( 7 ) ,NHEAD ( 20 ) , LABEL , ITAG , JTAG 

15020C***  THIS  SUBROUTINE  WRITES  ERROR  MESSAGES 

15030  NERORS=NERORS+l 

15040  GOTO  (10,20,30) ,NGOOF 

15050  10  WRITE  (LPRT,1001)  NCARD, LABEL 

15060  GO  TO  999 

15070  20  WRITE  (LPRT,1002)  NCARD , LABEL , ITAG , JTAG , ITAG 

15080  GO  TO  999 

15090  30  WRITE  (LPRT,1003)  NCARD , LABEL , ITAG , JTAG , LABEL 

15100  GO  TO  999 

15110  999  RETURN 

15120  1001  FORMAT  ('  IN  CARD  NO. ',15,'  OF  CURRENT  DATA  FILE,  THE  LABEL 
15130  &A4 , ' CAN  NOT  BE  INTERPRETED') 

15140  1002  FORMAT  ('  IN  CARD  NO.', 15,’  OF  CURRENT  DATA  FILE,  LABELLED 
15150  &A4 , IX , 11 , 13/ 1 THE  VALUE  OF  ITAG,  ',11, 

15160  &’  IS  NOT  IN  THE  PERMISSABLE  RANGE’) 

15170  1003  FORMAT  ('  CARD  NO.', 15,'  OF  CURRENT  DATA  FILE,  LABELLED 
15180  &A4 , IX , II , 13/'  IS  INCONSISTENT  WITH  A PREVIOUS  ', 

15190  &A4 , ' CARD') 

15200  9999  RETURN 
15210  END 

15220  SUBROUTINE  INTEG  (NEQ , X , Y , ISTAR, STPSZ ) 

15230C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

15240  COMMON/CMPRNT/IN ,LPRT,MONITR 

15250  COMMON/CMNTIG/FIRSTP, STEP, EPS, AB,NCUTS 

15260  INTEGER  NEQ,  NCUTS 

15270C***  ***  REAL*8  X , STEP , Y ( 6 ) , EPS , AB 

15280  LOGICAL  STPSZ 

15290C***  ***  REAL*8  HC/0 . 0D0/ , FINAL, H2 , H3 , H6 , H8 , ERR, TEST , T, H , EPSL , TEMPO 

15300  DIMENSION  Y (6 ) , Yl (6 ) , Y2 (6 ) , F0 (6 ) , FI ( 6 ) , F2  ( 6 ) 

15310  INTEGER  CUT 

15320  LOGICAL  DBL 

15330  DATA  HC/0.0/ 

15340C***  ***  REAL*8  Yl(6)  ,Y2(6)  ,F0(6)  , F 1 ( 6 ) ,F2(6) 

15350  50  FORMAT  ('  THE  STEPSIZE  IS  NOW' , 1PD15 . 6 , ' AT  TAU  =',D15.6) 

15360  60  FORMAT  ('  THE  STEPSIZE  HAS  BEEN  HALVED  ’,13,'  TIMES') 

15370  IF (NEQ .NE . 0)  GO  TO  10 

15380  HC  = FIRSTP 

15390  RETURN 

15400  10  IF (STEP. EQ.0)  RETURN 

15410  IF(HC.EQ.O)  HC  = STEP 

15420  FINAL  = X+STEP 

15430  H = STEP 

15440  EPSL  = EPS 

15450  IF (EPS . EQ. 0 .OR.ABS(H)  .LE.ABS(HC))  GO  TO  15 

15460  IF (H*HC . LE . 0)  HC  = -HC 

15470  H = HC 

15480  15  T = X+H 

15490  CUT  = NCUTS 
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15500  X = FINAL 

15510  H2  = H/2 . 

15520  H3  = H/3 . 

15530  H6  = H/6 . 

15540  H8  = H/8  . 

15550  20  I F ( H . GT . 0 .AND.  T.GT. FINAL  .OR.  H . LT . 0 . AND. T. LT . FINAL)  GOTO  40 

15560  21  CALL  FORCE (T-H ,Y , F0) 

15570  DO  22  I = 1 , NEQ 

15580  22  Y1 ( I ) = F0 ( I ) *H3+Y ( I ) 

15590  CALL  FORCE (T-2 . *H3 ,Yl , Fl) 

15600  DO  23  I = 1 , NEQ 

15610  23  Y1 ( I ) = (F0 (I)+Fl (I) ) *H6+Y ( I ) 

15620  CALL  FORCE (T-2 . *H3 ,Yl , Fl ) 

15630  DO  24  I = 1 ,NEQ 

15640  24  Yl ( I ) = (F1(I)*3.+F0(I))*H8+Y(I) 

15650  CALL  FORCE (T-H2,Yl ,F2) 

15660  DO  25  I = 1 ,NEQ 

15670  25  Yl ( I ) = (F2(I)*4.-F1(I)*3. +F0 ( I ) ) *H2  +Y(I) 

15680  CALL  FORCE (T,Y1 ,Fl) 

15690  DO  26  I = 1 ,NEQ 

15700  26  Y2 ( I ) = (F2 (I) *4 ,+Fl ( I) +F0 (I) ) *H6  +Y(I) 

15710  IF(EPSL.EQ.O)  GO  TO  38 

15720  DBL  = .TRUE. 

15730  DO  35  I = 1,NEQ 

15740  ERR  =ABS ( Yl ( I ) -Y2 ( I ) ) *0 . 2 

15750  TEST  =ABS (Yl (I) ) *EPSL 

15760  I F ( ERR . LE . TEST  .OR.  ERR. LT. AB)  GO  TO  34 

15770  H = H2 

15780  T = T-H2 

15790  IF  ( . NOT . STPSZ)  GO  TO  30 

15800  TEMPO  = T-H 2 

15810  WRITE  (LPRT , 50 ) H,  TEMPO 

15820  30  CUT  = CUT  - 1 

15830  IF  (CUT  .GE.  0)  GO  TO  31 

15840  X = T - H2 

15850  WRITE  (LPRT, 60)  NCUTS 

15860  RETURN 

15870  31  IF (T+H .NE . T)  GO  TO  33 

15880  X = T 

15890  RETURN 

15900  33  H2  = H/2. 

15910  H3  = H/3. 

15920  H6=  H/6. 

15930  H8  = H/8. 

15940  GO  TO  21 

15950  34  IF ( 64 . 0*ERR. GT . TEST)  DBL  = .FALSE. 

15960  35  CONTINUE 

15970  IF (.NOT. DBL)  GO  TO  38 

15980  H 2 = H 

15990  H = 2 . *H 
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16000 

IF  (STPSZ)  WRITE  (LPRT,50)  H , T 

16010 

H3  = H/3. 

16020 

H6  = H/6 . 

16030 

H8  = H/8 . 

16040 

CUT  = NCUTS 

16050 

38 

DO  39  I = 1 ,NEQ 

16060 

39 

Y(I)  = Y2  ( I ) 

16070 

T = T+H 

16080 

GO  TO  20 

16090 

40 

IF(EPSL.EQ.O)  RETURN 

16100 

HC  = H 

16110 

H = FINAL- (T-H) 

16120 

I F ( ABS ( H ) .LE.ABS (FINAL) *9 . 536744D-7 ) RETURN 

16130 

T=  FINAL 

16140 

EPSL  = 0 

16150 

H2  = H/2. 

16160 

H3  = H/3. 

16170 

H6  = H/6. 

16180 

H8  = H/8. 

16190 

GO  TO  20 

16200 

END 

16210 

SUBROUTINE  FORCE  ( TM , Y , DY ) 

16220C 

**  * 

***IMPLICIT  REAL*8  (A-H,0-Z) 

16230 

COMMON/CMPRNT/IN ,LPRT,MONITR 

16240 

COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

16250 

COMMON/CM LONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 

16260 

COMMON/CMCNST/RHO , GRAV , PI 

16270 

COMMON/CMDISP/TOWDSP , BOTDSP ,GYRAD , TMASS , ZNERTA , TLEN , TLEN2 , TLEN3 

16280 

COMMON/CMINIT/SPDIN ,CURENT , GAMIN , HEADIN ,CDOTIN , DRADIN , INSEG 

16290 

COMMON/CMSTAT/RPM (3) , SPEED , DELTA ( 2 , 3) , GAMMA , HEAD , CDOT , DRAD , 

16300 

&RPMHI (3) , ISEG 

16310 

COMMON /CMSI ZE/CGAFT , CGFWD , EFLEN , EFDRF , EFBEAM , TONS 

16320 

COMMON/CMCLER/CLIN , XLIN , CLOUT , XLOUT , XSIGN , CLBOW , BOWCL 

16330 

COMMON/ CMCOEF/A (10) ,B(10) ,AA(10) ,BB(10) ,COEF(5) 

16340 

COMMON/CMC LOG/NPOLY , NSTRIP , KORPOL , KORSTR 

16350 

COMMON/CMSCR1 /ALRCOM 

16360 

DIMENSION  Y (6 ) ,DY(6) 

16370 

LOGICAL  JUMP 

16380 

EXTERNAL  SIGN 

16390C 

*** 

THIS  SUBROUTINE  COMPUTES  HYDRODYNAMIC  FORCES 

16400 

DY ( 2 ) =Y ( 1 ) 

16410 

DY (4 ) =Y ( 3 ) 

16420 

DY  ( 6 ) =Y ( 5 ) 

16430 

CALL  CURRNT ( Y ( 2 ) ,Y(4)  ,Y(6)  ,CTAN ,CRAD,CROT) 

16440 

CY=Y ( 5 ) -CRAD 

16450 

CX=Y (6 ) *Y (1 ) 

16460 

SPD=SQRT(CY*CY+(CX*XSIGN+CTAN) **2) 

16470 

APGAMA=ATAN2 (CY, (CX+CTAN*XSIGN) ) 

16480 

IF  (XSIGN  .LT.  0.0)  APGAMA=APGAMA+PI 

16490 

APBETA=Y ( 4 ) -APGAMA 

16500  VADV=SPD*COS (APBETA) 

16510  ANGRAT=Y ( 3 ) -Y ( 1 ) -CROT 

16520  CALL  WINDFO  (Y ( 2) , Y ( 4 ) ,FWLON , FWLAT , FWROT) 

16530  CALL  THRUST  (VADV , PUSH , TORQ) 

16540  CALL  RESIST  (VADV, DRAG) 

16550  CALL  RUDFOR  (VADV , PROMOM , DELTA , FRUD , DRUD) 

16560  CALL  BOWTHR  (VADV , BTPUSH , BTSWAY , BTMOM) 

16570  JUMP=. FALSE. 

16580  CALL  HYDRO  (ANGRAT,SPD, APBETA, SMOM, FSWAY) 

16590  FPROP=PUSH-DRAG+BTPUSH+FWLON 

16600  10  CONTINUE 

16610  FSIDE=FRUD+FSWAY+BTSWAY+FWLAT 

16620  QALPH=FPROP*COS (Y (4) ) -FSIDE*SIN (Y (4) ) 

16630  QALPH=QALPH*Y (6) 

16640  QRAD=FPROP*SIN (Y (4) ) +FSIDE*COS (Y (4) ) 

16650  QHEAD=SMOM-FRUD* ( CGAFT+BOTLEN ) +BTMOM+PROMOM+ FWROT 

16660  RRM=TMASS*Y (6) *Y (6) 

16670  QKEEP=QALPH/RRM+QHEAD/RRM 

16680  DY(1)=QKEEP-2.0*Y(5) *Y(1)/Y(6) 

16690  DY (3 ) =QKEEP+QHEAD/ZNERTA-2 . 0*Y (5) *Y (1)  /Y (6) 

16700  DY ( 5 ) =QRAD/TMASS+Y ( 6 ) * Y ( 1 ) *Y ( 1 ) 

16710  IF  (JUMP)  GO  TO  100 

16720  ACC=SQRT(DY (1) *DY (1) *Y (6) *Y (6) +DY (5) *DY (5) ) 

16730  GA=ATAN2 (DY ( 5) ,Y (6) *DY (1) ) 

16740  BA=Y ( 4 ) -GA 

16750  VDOT=-ACC*SIN (BA) 

16760  FSWAY=FSWAY+B(10) *TLEN*TLEN*TLEN*VDOT*RHO/2 . 0 

16770  SMOM=SMOM+A ( 9 ) *TLEN*TLEN*TLEN*TLEN* ( DY ( 3 ) -DY ( 1 ) ) *RHO/2.0 

16780  JUMP= . TRUE . 

16790  GO  TO  10 

1G800  100  CONTINUE 

16810C***  WRITE  (MONITR, 1)  DRAG , PUSH , FPROP , FSWAY , FRUD , FSIDE 
16820C***  WRITE  (MONITR, 1)  QALPH ,QRAD,QHEAD,SPD, APGAMA 
16830C***  WRITE  (MONITR, 1)  (DY ( I ) , 1=1 , 6) 

16840C***  1 FORMAT  {’  CHECK 1 , 6E15 . 5) 

16850  RETURN 

16860  END 

16870  SUBROUTINE  PLOTER (RADIUS , PSI , ALPHA , I ) 

16880C***  ***IMPLICIT  REAL* 8 (A-H,0-Z) 

16890  COMMON/CMPRNT/IN , LPRT, MONITR 

16900  COMMON/CMDATA/TEMP ( 7 ) ,NHEAD (20 ) , LABEL , ITAG ,JTAG 

16910  COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

16920  COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 

16930  COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 

16940  COMMON/CMSTAT/RPM ( 3 ) , SPEED , DELTA ( 2 , 3)  , GAMMA, HEAD, CDOT,DRAD, 

16950  SRPMHI (3) , ISEG 

16960  COMMON/CMCNST/RHO ,GRAV , PI 

16970  COMMON/CMSIZE/CGAFT,CGFWD, EFLEN , EFDRF , EFBEAM ,TONS 

16980  COMMON/CMCLER/CLIN ,XLIN , CLOUT , XLOUT ,XSIGN ,CLBOW ,BOWCL 

16990  COMMON/CMSYMB/ROW ( 80 ) ,SYM(4) ,BL ,NMAX 
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17000 

17010 

17020 

17030 

17040 

17050 

17060 

17070 

17080 

17090 

17100 

17110 

17120 

17130 

17140 

17150 

17160 

17170 

17180 

17190 

17200 

17210 

17220 

17230 

17240 

17250 

17260 

17270 

17280 

17290 

17300 

17310 

17320 

17330 

17340 

17350 

17360 

17370 

17380 

17390 

17400 

17410 

17420 

17430 

17440 

17450 

17460 

17470 

17480 

17490 


COMMON/SCR2/NOUTS (8) ,POUTS(8) ,NSTRS,KNT 
DATASYM , BL/1HI , 1HO , 1HB , 1HT, 1H  / 

DATANOUTS/5HSPEED , 5HRPM-P , 5HRPM-S , 5H  -CL , 5HV-RAD , 5HV-YAW , 
& 5HS-RUD , 5HF-RUD/ 

DIST (X , Y) =SQRT(X*X+Y*Y+2.0*X*Y*SIN (PS I) ) 

RPRIM(Z) =RADIUS+Z/COS(PSI) 

XLEN (V , W) =V-W*SIN (PSI) /COS (PSI) 

EXES (T , U) =XSIGN* ( RADIUS ‘ALPHA +T*COS (PSI ) -U) 

NSTRS=NSTRS+1 

I F ( (NSTRS+1 ) /2*2 . EQ . (NSTRS+1 ) ) WRITE ( MONITR , 96 ) NHEAD 

96  FORMAT ('l' , 20A4) 

XSIGN=SANG ( 1 , I)/ABS (SANG (1 ,1) ) 

TWID=TOWWID/2 . 0*XSIGN 
BWID=BOTWID/2 . 0*XSIGN 
PCRT=ALPHA/SANG (1 , I) *5730 . 

WRITE (MONITR, 97 ) PCRT, I 

97  FORMAT ( //F5. 1, ’ PERCENT  THRU  SEGMENT', 13) 

XINC= (TOWLEN+BOTLEN) /7 . 

POUTS ( 1 ) =SPE ED/1 . 47 
POUTS ( 2 ) =RPM ( 1 ) 

POUTS ( 3 ) =RPM ( 2 ) 

POUTS ( 4 ) =RPM ( 3 ) 

POUTS (5) =SPEED*SIN (GAMMA) /I . 47 
POUTS ( 6 ) =CDOT*57 . 3 
POUTS ( 7 ) =DELTA ( 1 , 1) *57 . 3 
CLIN=9999 . E+10 
POUTS ( 8 ) =DELTA (2,1) *57.3 
CLOUT=9999 . E+10 
DO102JP=l ,25 
NMAX=0 
DO10K=l , 80 
10  ROW (K) =BL 

XPOS=XINC* FLOAT (20-JP) 

CALL  BARDIS(ALPHA, RADIUS, PSI, XPOS, DPERI ,DPERO, 

&PERRAD , PERTAN , ISSEGO , ISEGO) 

IF (DPERI .GT. CLIN) GOTO 12 

CLIN=DPERI 

XLIN=XPOS 

12  IF (DPERO.GT . CLOUT) GOTO) 4 
CLOUT=DPERO 
XLOUT=XPOS 
14  CONTINUE 

CALL  SCALE (-DPERI, 1,1) 

CALL  SCALE (DPERO ,1,2) 

IF (XPOS .GT.CGFWD) GOTOIOO 

I F ( XPOS . LT . -CGAFT-BOTLEN ) GOTOl 0 0 

WWID=TWID 

NSYM=3 

IF (XPOS . GT . -CGAFT) GOTO20 
WW,ID=BWID 
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1 7500 

17510  20 

17520 

17530 

17540 

17550 

17560  100 

17570 

17580 

17590  98 

17600 

17610  101 

17620  102 

17630  99 

17640 

17650 

17660 

17670C*** 

17680 

17690 

177  00 

1 7710 

1 .’72  J 

l^DO 

17  40 

17730 

177COC*** 

17770 

17780 

17790 

17800 

17  810 

1 7820 

17830 

17840  10 

17850 

17860 

17870C*** 

17880 

17890 

17900 

17910 

17920 

17930 

17940 

17950 

17960C*** 

17970C*** 

17980C*** 

17990 


NSYM=4 

CONTINUE 

CL=-WWID/COS (PSI) 

CALL  SCALE (CL , I ,NSYM) 

CL=WWID/COS (PSI) 

CALL  SCALE(CL,I,NSYM) 

JK=JP-14 

I F ( JK . LT . 1 . OR . JK . GT . 8 ) GOTOl 0 1 

WRITE ( MONITR ,98) NOUTS ( JK ) , POUTS (JK) , (ROW(JJ) ,JJ=11,NMAX) 

FORMAT ( IX, A5,F5.1,70A1) 

GOTOl 02 

WRITE (MONITR, 99 ) (ROW(JJ) ,JJ=11,NMAX) 

CONTINUE 
FORMAT ( 1 IX , 70A1 ) 

RETURN 

END 

SUBROUTINE  SCALE ( CL , I , ICHAR) 

***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON /CMPRNT/IN , LPRT , MONITR 
COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

COMMON /CM LONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 
COMMON/CMROUT/SDAT (11,10,2) ,CUR ( 11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 
COMMON/CMCNST/RHO , GRAV , PI 

COMMON/CMS I Z E/CG AFT , CGFWD , EFLEN , EFDRF , EFBE AM , TONS 
COMMON/CMCLER/CLIN , XLIN , CLOUT , XLOUT , XS IGN , CLBOW , BOWCL 
COMMON/CMSYMB/ROW ( 80 ) ,SYM(4) , BL,NMAX 
TEMPORARY  FIX  FOR  SEG  WIDTH 
SEGW=SDAT (1,1,2)  -SDAT (1,1,1) 

CL=XSIGN*CL 

N = IFIX (CL/SEGW*35 . ) +45 
I F (N . LT . 11 ) GOTOl 0 
I F (N .GT. 80 ) GOTOl 0 
ROW (N)=SYM( ICHAR) 

I F ( N . GT . NMAX ) NMAX=N 

RETURN 

END 

SUBROUTINE  THRUST  (VADV , PUSH , TORQ) 

♦♦♦IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT , MONITR 

COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD ,NPROP , NBLAD 
COMMON/CMCNST/RHO ,GRAV , PI 

COMMON/CMDISP/TOWDSP ,BOTDSP ,GYRAD ,TMASS , ZNERTA ,TLEN , TLEN2 , TLEN3 
COMMON/CMSTAT/RPM ( 3 ) ,SPEED , DELTA ( 2 , 3) , GAMMA , HEAD ,CDOT , DRAD , 

&RPMHI (3) , ISEG 
COMMON/SCR4/PUSHEA ( 3) 

LOGICAL  LIMIT, NIMIT 

THIS  SUBROUTINE  COMPUTES  THE  THRUST  AND  TORQUE  AT  THE  PRESEN 
PROPELLER  RPM  AND  ADJUSTS  THE  RPM  TO  THE  MAXIMUM  PERMITTED  B 
EITHER  COMMAND  OR  HORSEPOWER  LIMITATIONS 
PUSH=0 . 0 
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18000  DO  200  IP=1 ,NPROP 

18010  RPM (IP) =RPMHI (IP) 

18020  LIMIT=. FALSE. 

18030  NIMIT=. FALSE. 

18040  LOOP=0 

18050  XSIN=SIGN (1.0, RPMHI ( IP) ) 

18060  5 CONTINUE 

18070  RPS=RPM ( IP) /60 . 

18080  LOOP=LOOP+l 

18090  IF  (LOOP  .GT.  100)  GO  TO  50 

18100  IF  (NBLAD  .LT.  0)  GO  TO  20 

18110  CALL  QUAD4 (VADV , RPS , PUSHEA ( IP) ,TORQ) 

18120  GOTO  30 

18130  20  CONTINUE 

18140  C J= ( 1 . 0-WFRAC) *VADV/RPS/DPROP 

18150  CALL  WAGN4  (CJ,TKT,QKQ) 

18160  RN2D4=RHO*RPS*RPS*DPROP**4 

18170  PUSHEA (IP) =TKT*RN2D4* (1.0-TDDUC) 

18180  TORQ=QKQ*RN2D4*DPROP 

18190  30  CONTINUE 

18200  HP=PI*ABS(RPS*TORQ)/275.0 

18210  IF  (HP  .GT.  SHP  .OR.  ABS(RPM(IP))  .GT.  ABS (RPMHI ( IP) ) ) GO  TO  40 

18220  NIMIT= . TRUE . 

18230  IF  (LIMIT)  GO  TO  50 

18240  RPM (IP) =RPM (IP) +5 . *XSIN 

18250  IF  (ABS (RPM(IP) ) .LT.  ABS (RPMHI ( IP) ) +3 . ) GO  TO  5 

18260  RPM (IP) =RPM (IP) -5 . *XSIN 

18270  GO  TO  50 

18280  40  CONTINUE 

18290  IF (NIMIT) GO  TO  42 

18300  GO  TO  46 

18310  42  IF  (.NOT.  LIMIT)  GO  TO  45 

18320  GO  TO  46 

18330  45  NIMIT=. FALSE. 

18340  46  LIMIT= . TRUE . 

18350  IF  (NIMIT)  GO  TO  50 

18360  RPM ( IP) =RPM (IP)-5.*XSIN 

18370  GOTO  5 

18380  50  CONTINUE 

18390C***  WRITE  (LPRT, 100)  VADV, RPS , DPROP ,CJ 

18400C***  WRITE  (LPRT,100)  RPM ( IP) , RN2D4 , TKT ,QKQ 

18410C***  WRITE  (LPRT, 100)  PUSH , TORQ , HP , SHP 

18420C***  100  FORMAT  ('  PROP' ,4E15 . 5) 

18430  200  PUSH=PUSH+PUSHEA (IP) 

18440  RETURN 

18450  END 

18460  SUBROUTINE  RESIST  (VADV, DRAG) 

18470C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

18480  COMMON/CM  PRNT/IN , LPRT , MON I TR 

18490  COMMON/CMS I Z E/CGAFT , CGFWD , EFLEN , EFDRF , EFBEAM , TONS 
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18500 

18510C*** 

18520C*** 

18530C*** 

18540 

18550 

18560 

18570 

18580 

18590 

18600  10 
18610 
18620 

18630  20 

18640 
18650 
18660 
18670 
18680 
18690 
18700C*** 
18710 
18720 
18730C*** 
18740 
18750 
18760 
18770 
18780 
18790 
18800 
18810 
18820 
18830 
18840 
18850 
18860 
18870 
18880 
18890 
18900 
18910 
18920 
18930 
18940C*** 
18950 
18960 
18970 
18980 
18990 


EXTERNAL  SIGN 

THIS  SUBROUTINE  COMPUTES  DEEP  WATER,  UNRESTRICTED  CHANNEL  RE 
RESISTANCE 

FOR  THE  BARGE  AND  TOWBOAT 
EL2=SQRT (EFLEN) 

IF  ( EL 2 .GT.  20.0)  GO  TO  10 
CDRAG=1.723578*EL2 
CDRAG=CDRAG+5 . 520842E-3*EL2** 3 
CDRAG=CDRAG-2 . 542E-6*EL2**5 
GO  TO  20 
CONTINUE 

CDRAG=2.63833*EL2 
CDRAG=CDRAG+2 . 21667 E- 3* EL 2**3 
CONTINUE 

DRAG=184 .0* (TONS/CDRAG) **2 . 86 
DRAG=DRAG* (ABS (VADV) /EFBEAM) **1.86 
DRAG=DRAG*SIGN (1.0 ,VADV) /EFDRF**2 . 49 


RETURN 

END 

SUBROUTINE  WAGN4  (CJ , TKT,QKQ) 

***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN ,LPRT,MONITR 

COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD , NPROP , NBLAD 
COMPUTE  THRUST  COEFFICIENT 
TKT=-0 .719975E-2 
TKT=TKT-0.790916E-1*ARAT 
TKT=TKT-0. 179541  *CJ 
TKT=TKT-0.625748E-1*ARAT*CJ 
TKT=TKT-0. 311639  *CJ*CJ 
TKT=TKT+0. 143160  *ARAT*ARAT*CJ*CJ*CJ 

TKT=TKT+0 . 531326  *PD 
TKT=TKT-0. 114389  *ARAT*PD*CJ 
TKT=TKT+0 . 625376E-1*PD*CJ*CJ 
TKT=TKT+0. 125537  *PD*CJ*CJ*CJ 
TKT=TKT-0 . 52382 1E-1*ARAT*PD*CJ*CJ*CJ 
TKT=TKT-0 .207108  *PD*PD 
TKT=TKT+0. 270781  *ARAT*PD*PD 
TKT=TKT+0. 134182  *PD*PD*CJ 
TKT=TKT-0. 121086  *ARAT*PD*PD*CJ 
TKT=TKT-0 .189764E-1 * ARAT* ARAT* ARAT* PD* PD* CJ 
TKT=TKT-0 .439535E-1*ARAT*ARAT*ARAT*PD*PD*CJ*CJ 
TKT=TKT-0 . 624937E-1*PD*PD*CJ*CJ*CJ 
TKT=TKT-0 . 496939E-2* ARAT  * ARAT*  PD*  * 6 
TKT=TKT+0 . 115986E-1*ARAT*ARAT*PD**6*CJ 


COMPUTE  TORQUE  COEFFICIENT 
QKQ=0 . 964375E-2 
QKQ=QKQ-0 . 104 103E-1*ARAT 
QKQ=QKQ+0 . 512431E-2*ARAT*ARAT 
QKQ=QKQ+0 . 109936E-1 *ARAT*ARAT*ARAT 
QKQ=QKQ-0 . 453419 E-2*CJ 
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19000 

19010 

19020 

19030 

19040 

19050 

19060 

19070 

19080 

19090 

19100 

19110 

19120 

19130 

19140 

19150 

19160 

19170 

19180 

19190 

19200 

19210 

19220 

19230 

19240C*** 

19250 

19260 

19270 

19280 

19290 

19300 

19310 

19320 

19330 

19340 

19350 

19360 

19370 

19380 

19390 

19400 

19410 

19420 

19430 

19440 

19450 

19460 

19470 

19480 

19490 


QKQ=QKQ+0 . 216078E-1*ARAT*CJ 

QKQ=QKQ-0 . 507337E-1*CJ*CJ 

QKQ=QKQ+0 . 377970E-1*ARAT*CJ*CJ 

QKQ=QKQ-0.549486E-1*ARAT*ARAT*ARAT*CJ*CJ*CJ 

QKQ=QKQ-0 . 507  319E-1* ARAT*ARAT*PD 

QKQ=QKQ+0. 368649E-1*PD*CJ 

QKQ=QKQ-0. 106520  *ARAT*PD*CJ 

QKQ=QKQ+0 .465315E-1*ARAT*ARAT*ARAT*PD*CJ*CJ 

QKQ=QKQ+0 . 883010E-1*ARAT*ARAT*PD*CJ*CJ*CJ 

QKQ=QKQ+0 . 112619E-1*PD*PD 

QKQ=QKQ+0. 104825  *ARAT*PD*PD 

QKQ=QKQ-0 . 449154E-1*ARAT*PD*PD*CJ 

QKQ=QKQ+0 . 378780E-1*ARAT*ARAT*PD*PD*CJ 

QKQ=QKQ+0 . 177  304E-1*PD*PD*CJ*CJ 

QKQ=QKQ-0.164687E-1*ARAT*PD*PD*CJ*CJ 

QKQ=QKQ-0 . 344328E-1* ARAT*ARAT*PD*PD*CJ*CJ 

QKQ=QKQ-0 .249132E-1*ARAT*ARAT*ARAT*PD*PD*CJ*CJ 

QKQ=QKQ-0.233007E-1*ARAT*PD*PD*CJ*CJ*CJ 

QKQ=QKQ-0 .120209E-2*PD**6 

QKQ=QKQ-0 . 118997E-2*ARAT*ARAT*ARAT*PD**6 

QKQ=QKQ+0 . 458094E-2*ARAT*PD**6*CJ 

RETURN 

END 

SUBROUTINE  QUAD4 (VADV , RPS , PUSH1 , TORQ) 

★★★IMPLICIT  REAL* 8 (A-H ,0-Z) 

COMMON/CMPRNT/IN , LPRT , MONITR 

COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD , NPROP , NBLAD 

COMMON/CMCNST/RHO ,GRAV , PI 

DIMENSION  AT ( 21 ) ,BT(21) ,AQ(21) ,BQ(21) 

DATA  AT/0 .25350 E-1,0.17820E+0, 0.14674 E-l, 

&0.28054E-1,-. 16328 E-l, -.53041 E-l , 

&0 . 6060 5E-3 ,0.36823E-1  ,-.25429E-2, 
&-.17680E-l,0.27331E-2,0.21436E-l, 

&- . 24782E-2 ,0.1231 7E-2 ,0.50980E-2, 
&0.78076E-2,-.37816E-2,0.35353E-2, 

&0. 53014 E- 2,0. 21940E- 2,-. 28306 E-2/ 

DATA  BT/0.00000E+0 ,-.74777E+0 ,-.13822E-l , 

&0 . 10077E+0 , - . 11318E-1 ,0.47186 E-l , 

&0.10666E-1,-. 90 239 E-2, -.78452 E-2, 

&0. 23941 E-l,0.80787E-2,-. 14942 E-3 , 

&-. 3192 5E-2,0.92620E-2,0. 1552 7E-2, 

&-. 65683E-2, -.61655 E-3, 0.51033E-2, 

&-. 6026 3E-3,-. 82244 E-2,-. 63789E-3/ 

DATA  AQ/0. 2464 5E-1,0. 2671 8E+0 ,0.16056E-1 , 
&0.65822E-1,-.22497E-1,-. 78062 E-l , 

&0.24126E- 2, 0.61475 E-l, -.16065 E-l, 

&-.33  291E-1 ,0.1231 IE-1,0. 3112 3E-1 , 
&-.12559E-l,0.13948E-l,0.88397E-2  , 

&0. 50 358 E-3, -.79990 E-2, 0.13345E-1, 
&0.11928E-l,-.13556E-2,-.70825E-2/ 
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19500  DATA  BQ/0 . 00000E  + 0 , - . 11081E+1 , 0 . 15909E-2  , 

19510  &0. 164 5 5E+0 ,-. 20601E-1, 0.8534 3E-1 , 

19520  &0. 87856 E- 2, -.31327 E-l, -.96 650 E-2, 

19530  &0. 43190 E- 1,0. 12453E-1 ,0.95986E-3, 

19540  &-.79986E-2,0.15073E-l,0 .24595E-2 , 

19550  &-.16918E-l,0. 51603E-2 ,0.11504 E-l, 

19560  &-. 47976E-2,-. 14566 E-l, 0. 2328 OE-2/ 

19570  TKT=0 . 0 

19580  QKQ=0.0 

19590  WF=WFRAC 

19600  I F ( VADV . LT . 0 . ) WF=0 . 0 

19610  PFAC  = 0 . 7 *PI *RPS*DPROP 

19620  C JB=ATAN2 ( (1.0-WF) *VADV,PFAC) 

19630  DO  10  K=1 , 21 

19640  XK=FLOAT (K-l ) *CJB 

19650  TKT=TKT+AT (K) *COS (XK) +BT (K) *SIN (XK) 

19660  10  QKQ=QKQ+AQ (K) *COS (XK) +BQ(K) *SIN (XK) 

19670  PUSH1  = 0 . 125*RHO* (VADV*VADV+PFAC*PFAC)  *PI *DPROP*DPROP 

19680  TORQ=PUSH1*QKQ*DPROP/10 . 

19690  PUSH1=TKT*PUSH1 

19700  RETURN 

19710  END 

19720  SUBROUTINE  TIMCOR  ( Y , TSEG , KSTOP) 

19730C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

19740  COMMON/CMPRNT/IN , LPRT , MONITR 

19750  COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

19760  COMMON/CM INIT/SPDIN , CURENT , GAMIN , HEADIN , CDOTIN , DRADIN , INSEG 

19770  COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 

19780  COMMON/ CMCNST/RHO ,GRAV , PI 

19790  COMMON/ CMDISP/TOWDSP, BOTDSP , GYRAD , TMASS , ZNERTA , TLEN , TLEN 2 , TLEN  3 

19800  COMMON/CMSTAT/RPM { 3 ) , SPEED , DELTA ( 2 , 3) , GAMMA , HEAD ,CDOT ,DRAD , 

19810  &RPMHI ( 3) , ISEG 

19820  COMMON/CMSTER/BOWCLR, STRSLO , STRSLI , STRCOR , STRBK , RGAIN , RPMAX 

19830  DIMENSION  Y(6) 

19840C***  THIS  SUBROUTINE  COMPUTES  THE  ELAPSED  TIME  CORRECTION  REQUIRE 
19850C***  WHEN  THE  INTEGRATION  OVERSHOOTS  THE  END  OF  A SEGMENT 
19860C***  VARIABLES  DESCRIBING  THE  VELOCITY  AND  ORIENTATION  OF  THE  BAR 
19870C***  ARE  INITIALIZED  FOR  THE  START  OF  THE  NEXT  SEGMENT 
19880  IF  (KSTOP  .EQ.  4)  GO  TO  2 

19890  SPEED=SQRT (Y(1)*Y(1)*Y(6)*Y(6) +Y ( 5 ) *Y ( 5 ) ) 

19900  GAMMA=ATAN2 (Y(5) ,Y(6) *Y(1) ) 

19910  IF  ( SANG ( 1 , ISEG)  .LT.  0.0)  GAMMA=GAMMA+PI 

19920  SEGD= ( SDAT ( 1 , ISEG , 1 ) +SDAT ( 1 , ISEG , 2 ) ) /2 . 

19930  - EAD=Y ( 4 ) 

19940  DRAD=Y ( 6 ) -SEGD 

19950  CDOT=Y ( 3 ) 

19960  IF  (KSTOP  .EQ.  3)  GO  TO  999 

19970  TCOR=Y (6) * (ABS (Y (2) ) -ABS ( SANG ( 1 , ISEG) /57 . 3) ) /SPEED 

19980  HEAD=Y ( 4 ) -TCOR* Y ( 3 ) 

19990  DRAD=Y ( 6 ) -TCOR* Y ( 5 ) -SEGD 
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20000 

I=ISEG+1 

20010 

IF  (I  .GT.  NS EG)  GO  TO  30 

20020 

2 

CONTINUE 

20030 

CHSIGN=1 . 0 

20040 

J=ISEG 

20050 

5 

CONTINUE 

20060 

J=J  + 1 

20070 

IF  (J  .GT.  NS EG)  GO  TO  30 

20080 

SEGW=SDAT (1 , J, 2) -SDAT ( 1 , J , 1 ) 

20090 

IF  (SANG (1 , INS EG) *SANG ( 1 , J)  .LT.  0.0)  CHSIGN=-1.0 

20100 

HEAD=HEADIN*CHSIGN 

20110 

G AMMA=G AMIN * CHS IGN 

20120 

CDOT=CDOTIN*CHSIGN 

20130 

DRAD=DRADIN* ( BOWCLR*SEGW- (SEGW-TOWWID) /2 .0) 

20140 

IF  (SANG (1 , J)  .GT.  0.0)  GO  TO  30 

20150 

GO  TO  15 

20160 

10 

CONTINUE 

20170 

IF  (SANG (1,ISEG)*SANG (1,1)  .GT.  0.0)  GO  TO  30 

20180 

DRAD=-DRAD 

20190 

IF  (SANG (1,1)  .GT.  0.0)  GO  TO  20 

20200 

15 

CONTINUE 

20210 

HEAD=HEAD+PI 

20220 

GAMMA=GAMMA+PI 

20230 

GO  TO  30 

20240 

20 

CONTINUE 

20250 

HEAD=HEAD-PI 

20260 

GAMMA=GAMMA-PI 

20270 

30 

CONTINUE 

20280 

TSEG=TSEG-TCOR 

20290 

WRITE  (LPRT, 100)  SPEED, TCOR,CDOT, HEAD, GAMMA, DRAD 

20300 

100 

FORMAT  ('***S,T,C,H,G,D***' ,6F15.5) 

20310 

999 

RETURN 

20320 

END 

20330 

SUBROUTINE  RUDFOR ( VADV , PROMOM , DELTA , FRUD , DRUD) 

20340C* 

** 

***IMPLICIT  REAL* 8 (A-H,0-Z) 

20350 

COMMON/ CMPRNT/ IN , LPRT , MONITR 

20360 

COMMON/CM PRO P/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD ,NPROP , NBLAD 

20370 

COMMON/CMCNST/RHO ,GRAV , PI 

20380 

COMMON/CMRUDR/ARUD (2,3) , DELMAX , DLDTMX , OFSET 

20390 

COMMON/SCR4/PUSHEA ( 3 ) 

20400 

DIMENSIONDELTA (2,3) 

20410C*** 

THIS  SUBROUTINE  COMPUTES  RUDDER  LIFT  AND  DRAG  FORCES 

20420 

FRUD=0 . 0 

20430 

DRUD=0 . 0 

20440 

DO  100  IS=1 , 2 

20450 

DO  100  IT=1 ,NPROP 

20460 

CDRUD=0 . 0 

20470 

CDRUD=CDRUD+0 . 0*DELTA (IS , IT) 

20480 

CDRUD=CDRUD+0 . 0*DELTA ( IS , IT) *DELTA ( IS , IT) 

20490 

CDRUD=CDRUD+0 . 0*DELTA (IS , IT) *DELTA ( IS , IT) * DELTA ( IS , IT) 
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20500 
20510 
20520 
20530 
20540C 
20550 
20560 
20570 
20580 
20590 
20600 
20610 
20620 
20630 
20640C 
20650 
20660 
20670 
20680 
20690 
20700 
20710 
20720 
20730 
20740 
20750 
20760 
20770 
20780 
20790 
20800 
2081 0C 
20820C 
20830 
20840 
20850 
20860 
20870 
20880 
20890 
20900 
20910 
20920 
20930 
20940 
20950 
20960 
20970 
20980 
20990 


CLRUD=0 . 0 

CLRUD=CLRUD+2 . 0*DELTA ( IS , IT) 

CLRUD=CLRUD+0 . 0*DELTA ( IS , IT) *DELTA ( IS , IT) 

CLRUD=CLRUD+0  . Q*DELTA  ( IS,,  IT)  *DELTA  ( IS  , IT)  *DELTA  ( IS  , IT) 

***  TEMP  STMNT  TO  COMPUTE  AVG  PUSH 
APROP=PI*DPROP*DPROP/4 .0 

FACTOR=PUSHEA ( IT) /APROP+VADV* ABS (VADV) *RHO/2 . 0 
FRUD=CLRUD*ARUD ( IS , IT) * FACTOR+FRUD 
DRUD=CDRUD* ARUD ( IS , IT) *FACTOR+DRUD 
100  CONTINUE 

PROMOM= (PUSHEA ( 1 ) -PUSHEA ( 2 ) ) *OFSET 

RETURN 

END 

SUBROUTINE  DATOUT 
***  ***IMPLICIT  REAL* 8 (A-H ,0-Z) 

COMMON/CMPRNT/IN , LPRT ,MONITR 

COMMON/CMCNTR/NCARD ,NCASE ,NERORS , NPAGE , KPATH , KPRTNT , LINECT 
COMMON/CMDATA/TEMP (7 ) ,NHEAD (20 ) , LABEL , ITAG ,JTAG 
COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NBARL 
COMMON/CMCH AR/TOWDRF , BOTDRF , TOWBC , BOTBC 

COMMON/CMPROP/SHP , DPROP , PITCH , ARAT , WFRAC , TDDUC , PD ,NPROP , NBLAD 
COMMON/CMINIT/SPDIN ,CURENT , GAMIN , HEADIN ,CDOTIN , DRADIN , INSEG 
COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3 ) ,SANG(11,10) ,NANG(10) , NSEG 
COMMON/CMSTAT/RPM ( 3 ) , SPEED, DELTA (2 , 3) , GAMMA , HEAD , CDOT , DRAD , 

&RPMHI (3) , ISEG 

COMMON/CMNTIG/FIRSTP , STEP , EPS , AB , NCUTS 
COMMON/CMRUDR/ARUD (2,3) , DELMAX , DLDTMX , OFSET 

COMMON/CMSTER/BOWCLR, STRSLO , STRSLI , STRCOR , STRBK , RGAIN , RPMAX 
COMMON /CMBOWT/BTHRUS ( 7 ) , BTSPD ( 7 ) , BTMAX , BTPOS , BTGA IN , NBTSPD , IDBT 
LOGICAL  IDBT 

***  THIS  SUBROUTINE  PRINTS  OUT  THE  INPUT  DATA  DESCRIBING  THE  BAR 
***  AND  TOWBOAT  CHARACTERISTICS 
WRITE  (LPRT, 9000) 

NHE AD , NBARL , NBARW 

TOWLEN , BOTLEN , TOWWID , BOTWID , TOWDRF , BOTDRF 
TOWBC , BOTBC , CGTOW , CGBOT , TOWK , BOTK 
NPROP , SHP , RPMAX , OFSET 
NBLAD , DPROP , PITCH , ARAT , WFRAC , TDDUC 
ARUD (1,1) , ARUD (2,1) , DELMAX , DLDTMX 
BOWCLR , SPDIN , GAMIN , HEADIN , CDOTIN , DRADIN 
NCUTS , FIRSTP , STEP , EPS , AB 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


(LPRT, 1000) 
(LPRT, 2000) 
(LPRT,  3->00) 
(LPRT, 4000) 
(LPRT, 5000) 
(LPRT, 6000) 
(LPRT, 7000) 
(LPRT, 8000) 
(LPRT, 9000) 


1000  FORMAT  { 20X , 20A4////10X ,' CHARACTERISTICS  OF  TOWBOAT  AND', 13, 

S'  LONG  BY', 13,'  WIDE  BARGE  TOW' ///37X ,' TOW'  , 11X ,' BOAT V) 

2000  FORMAT  (15X, 'LENGTH  OVERALL ’ , F12 . 0 , FI 5 . 0/15X , ' WIDTH ' , F21 . 0 , FI 5 . 0/ 
&15X , 'DRAFT' , F23 . 2 , Fl 5 . 2/) 

3000  FORMAT  (15X, 'BLOCK  COEFFICIENT ', Fl2 . 3 , Fl 5 . 3/15X ,' L C G (FORWARD) ', 
&F12.1,F15.1/15X, 'GYRADIUS ' , F19 . 1 , F15 . 1///) 

4000  FORMAT  { 10X ,’ PROPULSION  AND  RUDDER  CHARACTERISTICS  AND  STEERING', 
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21000  & ' CRITERIA’ ///15X , 'NUMBER  OF  SHAFTS’, 114/ 

21010  &15X, 'HORSEPOWER  PER  SHAFT ', FI 1 . 0/1 5X MAXIMUM  RPM',F20.0/ 

21020  &15X, 'SHAFT  OFFSET  FROM  CL',F11.0/) 

21030  5000  FORMAT  (15X, 'BLADES  PER  PROPELLER ', I 10/1 5X ,' DIAMETER ', F24 . 1/ 

21040  &15X, 'PITCH' ,F27.1/15X, 'AREA  RATIO ’ , F24 . 3// 

21050  &15X, 'WAKE  FRACTION ', F22 . 4/1 5X ,' THRUST  DEDUCTION ’ ,F19 . 4/) 

21060  6000  FORMAT  (15X,'AREA  PER  STEERING  RUDDER F7 . 0/1 5X , 

21070  & ' AREA  PER  FLANKING  RUDDER'  , F7 . 0/15X MAXIMUM  RUDDER  ANGLE FI 4 . 3/ 

21080  &15X, 'MAXIMUM  RUDDER  RATE ' ,F15 . 3///) 

21090  7000  FORMAT  ( 10X ,' INITIAL  BARGE/TOWBOAT  VELOCITY  AND  ORIENTATION'/// 
21100  &15X , ' BOWCLR  SPEED  GAMMA  YAW  YAW  RADIAL'/ 

21110  S.42X, 'ANGLE  RATE  OFFSET'/ 

21120  &15X,F6.2,F9.2,1X,3F8.4,F8.3///) 

21130  8000  FORMAT  ( 10X INTEGRATION  CONTROL  PARAMETERS'/// 

21140  &15X , 'NCUTS  FIRST-STEP  STEP-SIZE  REL-ERROR  ABS-ERROR'/ 

•21150  &15X,I5,2F12.4,2F12.8) 

21160  9000  FORMAT  ('1') 

21170  RETURN 

21180  END 

21190  SUBROUTINE  TRPOUT 

21200C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

21210  COMMON/CMPRNT/IN , LPRT , MONITR 

21220  COMMON /CMCNTR/NCARD , NCASE ,NERORS ,NPAGE , KPATH , KPRINT > LINECT 

21230  COMMON/CMDAT A/TEMP (7) ,NHEAD ( 20 ) , LABEL , ITAG ,JTAG 

21240  COMMON/CM INI T/SPDIN ,CURENT ,GAMIN , HEADIN , CDOTIN , DRADIN , INSEG 

21250  COMMON/CMROUT/SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) , NSEG 

21260  COMMON/CMTIME/ACUMT ( 100 ) ,STIME(100) , NXTCRV, LSTCRV ,NXTREV , LSTREV 

21270  LOGICAL  NXTCRV, LSTCRV, NXTREV, LSTREV 

21280C***  THIS  SUBROUTINE  PRINTS  A SUMMARY  OF  ELAPSED  TIME  FOR  A TRIP 
21290  CALL  DATOUT 

21300  WRITE  (LPRT, 1000)  NHEAD, NCASE 

21310  DO  100  I=INSEG , NSEG 

21320  SEGD= (SDAT (1,1,1) +SDAT (1,1,2) ) / 2 - 

21330  SEGW=SDAT (1,1,2) -SDAT (1,1,1) 

21340  WRITE  (LPRT, 2000)  I , SEGD , SEGW, SANG ( 1 , I ) ,STIME ( I ) , ACUMT ( I ) 

21350  100  CONTINUE 

21360  1000  FORMAT  (///20X , 20A4/////1 OX ,' SUMMARY  OF  ELAPSED  TIME  FOR  CASE', 
21370  &I3//10X , ' SEGMENT  RADIUS  WIDTH  ANGLE  ', 

21380  & ' TIME (SEG)  TIME ( ACUM) ' //) 

21390  2000  FORMAT  (10X , 15 , 2X , 2F1 0 . 0 , F10 . 4 , 2F12 . 2) 

21400  RETURN 

21410  END 

21420  SUBROUTINE  SETY  (Y, ALPHA) 

21430C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

21440  COMMON/CMPRNT/IN, LPRT, MONITR 

21450  COMMON/CMROUT/SDAT (11 ,10 ,2)  ,CUR (11 , 10 , 3)  ,SANG(11,10)  ,NANG(10)  , NSEG 

21460  COMMON/CMSTAT/RPM (3) , SPEED , DELTA ( 2 , 3) , GAMMA, HEAD, CDOT,DRAD, 

21470  S.RPMHI  (3)  , ISEG 

21480  DIMENSION  Y(6) 

21490  SEGD= ( SDAT ( 1, ISEG, 1)+SDAT(1, ISEG, 2) ) /2 . 
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21500  Y ( 6 ) =SEGD+DRAD 

21510  Y ( 2 ) = ALPHA 

21520  Y ( 3) =CDOT 

21530  Y (4) =HEAD 

21540  Y (5) =SPEED*SIN (GAMMA) 

21550  Y (1) =SPEED/Y (6) *COS (GAMMA) 

21560  RETURN 

2i 570  END 

21580  SUBROUTINE  SETIC  (Y ,TSEG, ALPHA, NTIMES ,NEMORY) 

21590C***  ***IMPLICIT  REAL*8  (A-H,0-Z) 

21600  COMMON/CMPRNT/IN ,LPRT,MONITR 

21610  COMMON/CMCNTR/NC ARD , NC ASE , NERORS , NPAGE , KP ATH , KPRINT , LINECT 

21620  COMMON/CM IN IT/SPDIN ,CURENT , GAMIN ,HEADIN ,CDOTIN , DRADIN , INS EG 

21630  COMMON/CMSTAT/RPM ( 3 ) , SPEED , DELTA (2,3) ,GAMMA,HEAD,CDOT , DRAD, 

21640  &RPMHI (3) ,ISEG 

21650  COMMON/CMBEGN/BEGRPM ( 3 ) ,BEGSPD , BEGDEL (2 , 3) ,BEGGAM , BEGHED , BEGCDT , 

21660  &BEGDRD 

21670  COMMON/ CMSTER/BOWCLR , STRSLO , STRSLI , STRCOR , STRBK , RGAIN , RPMAX 

21680  COMMON/CMS I Z E/CGAFT , CGPWD , EFLEN , EFDRF , EFBEAM , TONS 

21690  COMMON/CMC LER/CLIN , XLIN , CLOUT , XLOUT , SIGN , CLBOW , BOWCL 

21700  COMMON/CMTIME/ACUMT ( 100 ) ,STIME(100) ,NXTCRV,LSTCRV,NXTREV,LSTREV 

21710  COMMON/SCR2/NOUTS ( 8 ) ,POUTS(8) ,NSTRS ,KNT 

21720  COMMON/SCR3/HISTRY (10) 

21730  DIMENSION  Y(6) 

21740  LOGICAL  IMPROV,JUMP 

21750  LOGICAL  NXTCRV, LSTCRV ,NXTREV , LSTREV 

21760  MSIZE=6 

21770  CALL  TIMCOR  (Y,TSEG,3) 

21780  NPOINT=NTIMES- (NTIMES-1 ) /MSIZE*MSIZE 

21790  IF (DELTA (1,1) .GT . 1 . 57 ) GOTO40 

21800  IF (DELTA (1,1) .LT.-l . 57)G0T015 

21810  HISTRY (1) =TSEG 

21820  HISTRY ( 2 ) =ALPH A 

21830  HISTRY ( 3 ) =RPM (1 ) 

21840  HISTRY (4 ) =RPM (2) 

21850  HISTRY ( 5 ) =RPM ( 3 ) 

21860  HISTRY (6) =SPEED 

21870  HISTRY (7 )=GAMMA 

21880  HISTRY (8 )=HEAD 

21890  HISTRY (9) =CDOT 

21900  HISTRY (10) =DRAD 

21910  GO  TO  100 

21920  15  JUMP= . TRUE . 

21930  TSEG=0 . 0 

21940  ALPHA=0 . 0 

21950  SPEED=BEGSPD 

21960  GAMMA=BEGGAM 

21970  HEAD=BEGHED 

21980  CDOT=BEGCDT 

21990  DRAD=BEGDRD 
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22150 
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22280 

10 

22290 

22300 
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20 
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40 
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50 
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22400 

22410 

1000 

22420 

22430 

22440 

22450 

22460 

22470 

22480C*** 

22490 

GO  TO  50 
CONTINUE 
TSEG=HISTRY ( 1 ) 

ALPHA=HISTRY (2) 

RPM ( 1 ) =HISTRY ( 3 ) 

RPM ( 2 ) =H I STRY ( 4 ) 

RPM ( 3 ) =HISTRY (5) 

SPEED=HISTRY (6) 

GAMMA=HISTRY (7) 

HEAD=HISTRY (8) 

CDOT=HISTRY (9) 

DRAD=HISTRY ( 10 ) 

CONTINUE 

CALL  SETY  (Y, ALPHA) 

CALL  CHANCK  (Y (6 ) , Y (4 ) , Y ( 2 ) , ISEG) 

CONTINUE 

RETURN 

END 

SUBROUTINE  POLY 3 (XVAL, YVAL,COEF) 

♦♦♦IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/ CM PRNT/IN , LPRT , MONITR 

DIMENSION  X (6) ,XY(6) ,XVAL(3) ,YVAL(3) ,COEF(3) 

EQUIVALENCE  (X(l)  ,Xl)  , (X ( 2)  ,X2)  , (X ( 3)  ,X3) 

EQUIVALENCE  (X ( 4 ) ,X4)  , (X(5)  ,X5)  , (X(6)  ,X6) 

EQUIVALENCE  (XY(1)  ,X1Y)  , (XY(2)  ,X2Y)  , (XY(3)  ,X3Y) 

DO  10  1=1,6 

X(I)=0.0 

XY  (I ) =0 . 0 

CONTINUE 

DO  30  1=2,6 

DO  20  J=1 , 3 

X ( I ) =X ( I ) +XVAL ( J) * *1 

CONTINUE 

CONTINUE 

DO  50  1=1,3 

DO  40  J=1 ,3 

XY ( I ) =XY ( I ) +YVAL ( J) *XVAL (J) **I 


I 

i 


I 


CONTINUE 

CONTINUE 

D=X2*(X4*X6-X5*X5)-X3* (X3*X6-X4*X5) +X4* (X3*X5-X4*X4) 

IF  (D  .EQ.  0.0)  WRITE  (LPRT, 1000) 

FORMAT  ('DETERMINANT  IS  ZERO  IN  POLY 3 ' ) 

COEF (1)=(X1Y*(X4*X6-X5*X5) -X2Y* (X3*X6-X4*X5) +X3Y* (X3*X5-X4*X4) ) /D 
COEF (2) = (-X1Y* (X3*X6-X4*X5) +X2Y* (X2*X6-X4*X4) -X3Y* (X2*X5-X3*X4) ) /D 
COEF (3) = (X1Y* (X3*X5-X4*X4) -X2Y* (X2*X5-X3*X4) +X3Y* (X2*X4-X3*X3) ) /D 
RETURN 
END 

SUBROUTINE  SELECT  (PARAM , VALUE ,NEND , LIM) 

♦♦♦IMPLICIT  REAL* 8 (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT, MONITR 
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22500 

22510 

22520 
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22540 

22550 

22560 

22570 

22580 

22590 
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22600 

22610 

22620 

22630 

22640 

22650 

10 

22660 

22670 

22680 

22690 

20 

22700 

22710 

22720 

25 

22730 

22740 

22750 

22760 

22770 

30 

22780 

22790 

22800 

22810 

22820 

40 

22830 

22840 

22850 

22860 

22870 

22880 

50 

22890 

22900 

22910 

100 

22920 

22930 

22940 

22950 

1000 

22960 

2000 

22970 

22980 

22990 

DIMENSION  PARAM ( 7 ) , LIM ( 3) 

LIMIN=1 

LIMAX=NEND 

CHECK  *PARTIM ( I , J , K)  * FQR  ((-1)),  AND  ADJUST  LIMIN  AND  LIMAX 
PMIN=PARAM (LIMIN) 

PMAX=PARAM (LIMAX) 

IF  ( PMIN  .LE.  VALUE  .AND.  PMAX  .GT.  VALUE)  GO  TO  5 
WRITE  (LPRT, 1000)  VALUE , PMIN , PMAX 
GO  TO  100 
CONTINUE 

DO  20  I=LIMIN, LIMAX 
IBACK=LIMAX+1-I 

IF  (PARAM (I)  .GT.  VALUE)  GO  TO  10 
PLO=PARAM ( I ) 

MINN=I 

CONTINUE 

IF  (PARAM ( IBACK)  .LE.  VALUE)  GO  TO  20 
PHI=PARAM ( IBACK) 

MAXX=IBACK 

CONTINUE 

IF  ( MAXX  .EQ.  MINN+1 ) GO  TO  25 

WRITE  (LPRT, 2000)  MINN , PARAM ( MINN) , MAXX , PARAM (MAXX) 

CONTINUE 

IF  (MINN  .GT.  1)  GO  TO  30 

NXTRA=2 

MAXX=3 

GO  TO  100 

CONTINUE 

IF  (MAXX  .LT.  LIMAX)  GO  TO  40 

NXTRA=LIMAX-1 

MINN=LIMAX-2 

GO  TO  100 

CONTINUE 

AVG= ( PARAM ( MINN )+PARAM( MAXX) ) /2 . 0 

IF  (AVG  .GT.  VALUE)  GO  TO  50 

NXTRA=MAXX 

MAXX=MAXX+1 

GO  TO  100 

CONTINUE 

NXTRA=MINN 

MINN=MINN-1 

CONTINUE 

LIM (1 ) =MINN 

LIM (2) =NXTRA 

LIM ( 3 ) =MAXX 

FORMAT  (’  VALUE' ,F15.5, ' IS  OUT  OF  RANGE ' , 2F15 . 5) 

FORMAT  ('  PARAMETERS  OUT  OF  SEQUENCE ' , 2 ( 1 5 , 2X , F15 . 5) ) 

RETURN 

END 

SUBROUTINE  HYDRO  ( ANGRAT, SPD , APBETA , SMOM , FSWAY) 


BARGE 


PAGE  27 


23000C*** 
23010 
23020 
23030 
23040 
23050 
23060 
23070 
23080 
23090 
23100 
23110 
23120 
23130 
23140 
23150 
23160 
23170 
23180 
23190 
23200 
23210 
23220C*** 
23230 
23240 
23250 
23260 
23270 
23280 
23290 
23300 
23310 
23320 
23330 
23340 
23350 
23360 
23370 
23380 
23390 
23400 
23410 
23420 
23430 
23440 
23450 
23460 
23470 
23480 
23490 


100 


***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT ,MONITR 

COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NBARL 
COMMON/CMCNST/RHO ,GRAV , PI 

COMMON/CMDISP/TOWDSP ,BOTDSP ,GYRAD ,TMASS , ZNERTA, TLEN , TLEN2 , TLEN 3 
COMMON/CMCOEF/A (10) ,B(10) ,AA(10) ,BB(10) ,COEF(5) 

COMMON/ CMCLOG/NPOLY , NSTRI P , KORPOL , KORSTR 

TLEN4=TLEN*TLEN3 

HRO=RHO/2 . 0 

V=-SPD*SIN ( APBETA) 

Tl=HRO*TLEN2*SPD*SPD 
T2=HRO*TLEN2*SPD*V 
T3=HRO*TLEN3*SPD* ANGRAT 
T5=HRO*TLEN3/SPD* ANGRAT*V*V 
T6=HRO*TLEN4/SPD*ANGRAT* ANGRAT* V 
T7=HRO*TLEN2/SPD*V*V*V 

SMOM=A ( 1 ) *T1+A ( 2 ) *T2+A ( 3 ) *T3+A (5) *T5+A(6) *T6+A ( 7 ) *T7 
SMOM=SMOM*TLEN 

FSWAY=B(1) *T1+B(2) *T2+B (3) *T3+B (5) *T5+B (6) *T6+B (7) *T7 

RETURN 

END 

SUBROUTINE  HYCOEF 
***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMPRNT/IN , LPRT , MONITR 
COMMON /CMWIDE/TOWWID , BOTWID , NBARW 

COMMON/CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NBARL 
COMMON/CMCHAR/TOWDRF , BOTDRF , TOWBC , BOTBC 
COMMON/CMCNST/RHO ,GRAV , PI 

COMMON/CMDISP/TOWDSP ,BOTDSP ,GYRAD ,TMASS , ZNERTA , TLEN , TLEN2 , TLEN3 

COMMON/CMS I Z E/CGAFT , CGFWD , EFLEN , EFDRF , EFBEAM , TONS 

COMMON/CMCOEF/A (10) ,B(10) ,AA(10) ,BB(10) ,COEF(5) 

COMMON/CMCLOG/NPOLY , NSTRI P , KORPOL , KORSTR 

LOGICAL  NPOLY , NSTRIP , KORPOL , KORSTR 

TLEN=TOWLEN+BOTLEN 

TLEN2=TLEN*TLEN 

TLEN3=TLEN*TLEN2 

TLEN4=TLEN*TLEN3 

TLEN5=TLEN*TLEN4 

IF  (NPOLY)  GO  TO  800 

IF  (NSTRIP)  GO  TO  100 

IF  (KORPOL)  GO  TO  900 

IF  (KORSTR)  GO  TO  200 

CDT=1 . 5 

CDB=0 . 5 

CAT=0 .45 

CAB=0 .45 

CAS=0 . 0 

GO  TO  300 

CONTINUE 

KORPOL® . FALSE . 
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23500  200 

23510 

23520 
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23540 

23550 

23560  300 

23570 

23580 
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23640 

23650 

23660 

23670 

23680 
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23700 

23710 

23720  800 

23730 

23740  900 

23750 
23760 
23770 
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23790  1000 

23800 
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23820 
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23840 

23850 

23860 

23870 

23880 

23890 

23900 

23910 

23920 

23930 

23940 

23950  5 

23960 

23970 

23980 

23990 


CONTINUE 
CDT=COEF (1) 

CDB=COEF (2) 

CAT=COEF ( 3 ) 

CAB=COEF ( 4 ) 

CAS=COEF ( 5 ) 

CONTINUE 

ST=CAT*TOWBC*TOWWID*TOWDRF 
SB=CAB*BOTBC*BOTWID*BOTDRF 
DB=CGAFT+BOTLEN/2 . 0 
DT=CGFWD-TOWLEN/2 . 0 
TT  =TOWDRF  * TOWLEN 
BO=BOTDRF*  BOTLEN 

CC  = (CGAFT*CGAFT*CGAFT+CGFWD*CGFWD*CGFWD) /3 . 0 
CBC= ( -CGAFT*CGAFT*CGAFT+ (CGAFT+BOTLEN) **3) /3.0 
B ( 2 ) =-  (CDT*TT+CDB*BO) /TLEN2 

A ( 3) =- (CDT*TOWDRF*CC+CDB*BOTDRF*CBC) /TLEN4 
A ( 2) = (CDB*BO*DB-CDT*TT*DT) /TLEN3 
B ( 3 ) = A ( 2 ) 

A ( 9 ) =- ( ST*CC+SB*CBC) /TLEN5 

B (10) =- ( ST*TOWLEN+SB*BOTLEN) /TLEN3 

GO  TO  1000 

CONTINUE 

KORSTR=. FALSE. 

CONTINUE 
DO  950  K=1 , 10 
A (K) =AA (K) 

B (K) =BB (K) 

CONTINUE 

CONTINUE 

RETURN 

END 

SUBROUTINE  BOWTHR  (VADV ,BTPUSH ,BTSWAY ,BTMOM) 

♦♦♦IMPLICIT  REAL*8  (A-H.O-Z) 

COMMON /CMPRNT/IN , LPRT , MONITR 

COMMON/CMS I Z E/CGAFT , CGFWD , EFLEN , EFDRF , EFBEAM , TONS 
COMMON/CMCLER/CLIN , XLIN , CLOUT , XLOUT , SIGN , CLBOW , BOWCL 
COMMON/CMBOWT/BTHRUS (7 ) ,BTSPD(7) ,BTMAX ,BTPOS ,BTGAIN , NBTSPD , IDBT 
LOGICAL  IDBT 

DIMENSION  BTCOEF ( 3 ) ,BTH(3) ,V(3) ,LIMB(3) 

IF  (IDBT)  GO  TO  5 
BTPUSH=0 .0 
BTSWAY=0 .0 
BTMOM=0 . 0 
GO  TO  1000 
CONTINUE 

CALL  SELECT  (BTSPD, VADV, NBTSPD, LIMB) 

DO  10  1=1,3 
V (I) = BTSPD (LIMB (I ) ) 

BTH ( I ) =BTMAX-BTHRUS (LIMB ( I ) ) 
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CONTINUE 

CALL  POLY3  (V , BTH , BTCOEF) 

BTFORC=0 . 0 
DO  20  K=1 , 3 

BTFORC=BTFORC+BTCOEF (K) *VADV**K 
CONTINUE 

BTFORC=BTMAX-BTFORC 

IF  ( BTFORC  .GT.  0.0)  GO  TO  30 

BTPUSH=BTFORC 

BTSWAY=0 .0 

BTMOM=0 . 0 

GO  TO  1000 

CONTINUE 

BTLIM=BTMAX*BTGAIN* ( BOWCL-CLBOW) *SIGN/BOWCL 
IF  ( BTLIM  .GT.  BTMAX)  BTLIM=BTMAX 
IF  (BTLIM  .LT.  -BTMAX)  BTLIM=-BTMAX 
IF  (ABS (BTLIM)  .LT.  BTFORC)  GO  TO  100 
BTPUSH=0 . 0 

BTSWAY=BTFORC*BTLIM/ABS (BTLIM) 

BTMOM=BTSWAY* (CGFWD+BTPOS) 

GO  TO  1000 

CONTINUE 

BTSWAY=BTLIM 

BTMOM=BTSWAY* (CGFWD+BTPOS) 

BTPUSH=SQRT (BTFORC* BTFORC-BTS WAY* BTS WAY) 

CONTINUE 

RETURN 

END 

SUBROUTINE  DECEL  ( FRMSPD , TOSPD , TDIF) 

***IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON/CMNTIG/FIRSTP , STEP , EPS , AB , NCUTS 

COMMON/CMDISP/TOWDSP ,BOTDSP ,GYRAD , TMASS , ZNERTA , TLEN , TLEN2 , TLEN3 

TDIF=0 .0 

DSLOW=0 . 0 

NTRY=0 

BADS=FRMSPD 

IF  (FRMSPD  .LE.  TOSPD)  GO  TO  999 
CONTINUE 

CALL  RESIST  (BADS, DRAG) 

NTRY=NTRY+1 

BADT=BADS-DRAG*FIRSTP/TMASS 
DSLOW=DSLOW+FIRSTP* (BADS+BADT) /2 . 0 
BADS=BADT 

IF  (BADS  .GT.  TOSPD  .AND.  NTRY  .LT.  200)  GO  TO  10 

TSLOW=FIRSTP*NTRY 

TD I F=TSLOW-DS LOW/ FRMSPD 

RETURN 

END 

SUBROUTINE  BARDIS (ALPHA , RADIN , PSI , XPOS , DPERI ,DPERO, 

&PERRAD , PERTAN , ISSEG0 , ISEG0 ) 
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24500 
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24580 
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40 
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80 
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COMMON/ CMSTAT/DUM ( 17 ) , ISEGIN 

COMMON/CMROUT/SDAT (11,10,2)  ,CUR ( 1 1 , 10 , 3 ) ,SANG(11,10)  , NANG (10)  , NSEG 
DIMENSION  DISTI (10) ,DISTO(10) 

RADIUS=RADIN 
ISEG=ISEGIN 
DISTI (l)=9.999E+3 
DISTO (l)=9.999E+3 
PI180=3. 1415927/180. 

RAD=ALPHA 

CALCULATE  BARGE  POSITION  IN  SEGMENT  "ISEGIN" 

XGBAR=RADIUS*COS (RAD) 

YGBAR=RADIUS*SIN (RAD) 

CALCULATE  POINT  POSITION  IN  SEGMENT  "ISEGIN" 

RAD= (1 . 5708-PSI+ALPHA) 

IF ( ABS (XPOS) . LE . 0 . 001 ) XPOS=0 . 001 *SIGN (1.0, XPOS) 

XP=XGBAR+XPOS*COS (RAD) 

YP=YGBAR+XPOS*SIN (RAD) 

CALCULATE  ANGLE  OF  POINT  WITHIN  SEGMENT 
ANGPT=ATAN2 (YP,XP)/PI180 
ANGSEG=SANG ( 1 , ISEG) 

CHECK  TO  SIGN  OF  SEGMENT  ANGLE 
IF ( ANGSEG)  10,20,30 
ANG1=A0N2PI (ANGSEG) 

ANG2=360 . 

ANGPT0=A0N2PI (ANGPT) 

GO  TO  40 
RETURN 
ANG1=0 . 0 

ANG2=A0N2PI (ANGSEG) 

ANGPT0= ANGPT 

CHECK  IF  POINT  IS  WITHIN  SEGMENT  "ISEGIN" 

ISADD=0 

IF ( ANGPT0-ANG1)  50,80,60 

ISADD=-1 

GO  TO  80 

IF ( ANGPT0-ANG2)  80,80,70 

ISADD=1 

CONTINUE 

IF  BARGE  POINT  IS  NOT  IN  "ISEGIN" 

CG  IN  NEW  SEGMENT  "ISEG+ISADD" 

IF ( ISADD)  90,100,90 

CALL  COORD (XP,YP, ISEG, ISADD, XPI ,YPI) 

CALL  COORD ( XGBAR , YGBAR , I S EG , I S ADD , XG , YG ) 

XP®XPI 
YP=YPI 

WRITE (5,982)XP,YP 

982  FORMAT ( ' XP=',F8.2,'  YP=’,F8.2) 

XGBAR=XG 
YGBAR=YG 
IS EG® ISEG+ISADD 


DEFINE  POINT  AND  BARGE 


BARGE 


PAGE  31 


25000 

25010 

25020 

25030 

25040C 

25050 

25060 

25070 

25080 

25090 

25100 

25110 

25120 

25130 

25140 

25150 

25160C 

25170 

25180 

25190 

25200 

25210 

25220 

25230 

25240 

25250 

25260 

25270 

25280 

25290 

25300 

25310 

25320 

25330 

25340 

25350 

25360 

25370 

25380 

25390 

25400 

25410 

25420 

25430 

25440 

25450 

25460 

25470 

25480 

25490C 


ANGPT=ATAN2 (YP , XP) /PI 180 
ANGPT0=A0N2PI (ANGPT) 

RADIUS=SQRT (XP*XP+YP*YP) 

100  CONTINUE 

FIND  SUB-SEGMENT  WITHIN  "ISEG"  THAT  POINT  IS  LOCATED 
ANG=0 . 0 

NSSEG=NANG (ISEG) +1 
SN=SIGN ( 1 . , SANG ( 1 , ISEG) ) 

ANGPT=ANGPT0+1 80 .* (SN-1 . ) 

DO  110  ISS=2,NSSEG 
ISSEG0=ISS 

PERTAN= ( ANGPT-ANG) /SANG ( ISS , ISEG) 

ANG=ANG+SANG ( ISS , ISEG) 

I F(SN* (ANGPT-ANG) ) 120,120,110 
110  CONTINUE 
120  CONTINUE 

CALCULATE  INSIDE  AND  OUTSIDE  PERPENDICULAR  DISTANCES 
I1=MAX0 (ISEG-1 ,1) 

I2=MIN0 (ISEG+1 ,NSEG) 

10=0 

11=0 

DO  190  IS  = 11,12 
SN=SIGN ( 1 . , SANG (1 , IS) ) 

ISADD=IS-ISEG 

CALL  COORD (XP,YP, ISEG, ISADD,X,Y) 

CALL  COORD ( XGBAR , YGBAR , ISEG , ISADD , XG , YG ) 

DY=XG-X 

DX=Y-YG 

NSSEG=NANG (IS) +1 
X2I  = SDAT ( 1 , IS , 1 ) 

Y2I  = 0.0 

X20  = SDAT (1, IS, 2) 

Y20  =0.0 
ANG=0 . 0 

DO  190  ISS  = 2 , NSSEG 

XII  = X2I 

Yll  = Y2I 

XIO  = X20 

YlO  = Y20 

ANG  = ANG  + SANG (ISS, IS) 

R2I  = SDAT (ISS , IS , 1 ) 

R20  = SDAT (ISS , IS , 2) 

RAD  = ANG*PI 180 
COR  = COS (RAD) 

SIR  = SIN (RAD) 

X2I  = R2I*C0R 
Y2I  * R2I*SIR 
X20  = R20*COR 
Y20  * R20*SIR 

INSIDE  INTERSECTION  CALCULATIONS 
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25500 

25510 

25520 

25530 

25540C 

25550 

25560 

25570 

25580 

25590 

25600 

25610 

25620 

25630 

25640 

25650 

25660 

25670 

25680 

25690 

25700 

25710 

25720 

25730 

25740 

25750 

25760 

25770 

25780 

25790 

25800 

25810 

25820C 

25830C 

25840 

25850C 

25860 

25870 

25880C 

25890 

25900 

25910 

25920 

25930 

25940 

25950 

25960 

25970 

25980 

25990C 


CALL  LIN I NT (Xll , Yll , X2I , Y2I , X , Y , DX , DY , DIST , ISTAT) 

IF(ISTAT)  140,140,130 
130  11=11+1 

DISTI (II) =SN*DIST 
OUTSIDE  INTERSECTION  CALCULATIONS 
140  CALL  LININT (XlO  , YIO ,X20 , Y20 ,X , Y , DX , DY , DIST , ISTAT) 

IF (ISTAT)  160,160,150 
150  IO  = IO  + 1 

DISTO ( IO) =-SN*DIST 
160  CONTINUE 

IF(ISADD)  190,170,190 
170  IF (ISS-ISSEG0)  190,180,190 
180  CONTINUE 

CALL  LININT (XI I ,Y1I,X2I,Y2I,X,Y,X,Y,DIN, ISTAT) 

CALL  LININT (XlO, YIO, X20,Y20,X,Y,X,Y,D0UT, ISTAT) 

PERRAD=DIN/ (DIN-DOUT) 

190  CONTINUE 
ISEG0=ISEG 
DPERI=DISTI (1) 

DPERO=DISTO ( 1 ) 

IO=MAX0 (10,1) 

II=MAX0 (II  ,1) 

IF(IO.EQ.l)GOTO  201 
DO  200  1=2,10 

200  I F ( ABS ( DI STO ( I ) ) . LT . ABS (DPERO) ) DPERO=DISTO (I ) 

201  IF (II .EQ.l)GOTO  203 
DO  202  1=2,11 

202  I F ( ABS (DISTI (I) ) . LT . ABS (DPERI) ) DPERI=DISTI (I) 

203  CONTINUE 
RETURN 
END 

SUBROUTINE  COORD (XP, YP, ISEG, ISADD,XPI ,YPI) 

SUBROUTINE  CHANGES  COORDINATE  SYSTEM  OF  POINT  (XP,YP)  IN  SEGMENT 
"ISEG"  TO  (XPI.YPI)  IS  SEGMENT  "ISEG+ISADD" 

COMMON/CMROUT/SDAT (11,10,2) ,CUR ( 11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 
CHECK  WHICH  SEGMENT  FOR  POINT 
PI180=3. 1415927/180. 

IF(ISADD)  300,200,100 

CALCULATE  POINT  COORD  IN  SEGMENT  FORWARD  OF  BARGE 
100  ALl=SANG ( 1 , ISEG) 

NNS=NANG ( ISEG) +1 
RAD1=AL1 

IF (ALl*SANG (1 , ISEG+1) ) 110,500,120 
110  DR=SDAT (NNS , ISEG , 1 ) +SDAT ( 1 , ISEG+1 ,2) 

ROT=180 . +AL1 
GO  TO  400 

120  DR=SDAT (NNS , ISEG , 1 ) -SDAT ( 1 , ISEG+1 , 1 ) 

ROT=ALl 
GO  TO  400 

LEAVE  POINT  COORDS  IN  PRESENT  SEGMENT 
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26000 

200 

XPI=XP 

26010 

YPI=YP 

26020 

RETURN 

26030C 

CALCULATE  POINT  COORDS  IN  SEGMENT 

26040 

300 

AL1=SANG ( 1 , ISEG-1 ) 

26050 

NNS=NANG ( ISEG-1 ) +1 

26060 

RAD1=0 . 0 

26070 

I F ( AL1*SANG (1 , IS EG) ) 310,500,320 

26080 

310 

DR=SDAT (NNS , ISEG-1 , 1 ) +SDAT ( 1 , ISEG , 

26090 

ROT=180 . -AL1 

26100 

GO  TO  400 

26110 

320 

DR=SDAT ( 1 , ISEG , 1 ) -SDAT (NNS , ISEG-1 , 

26120 

ROT=-ALl 

26130 

400 

RAD=RAD1*PI180 

26140 

DH=DR*COS (RAD) 

26150 

DK=DR*SIN (RAD) 

26160 

XPT=XP-DH 

26170 

YPT=YP-DK 

26180 

RAD=ROT*PI 180 

26190 

COR=COS (RAD) 

26200 

SIR=SIN (RAD) 

26210 

XPI=XPT*COR+YPT*SIR 

26220 

YPI=YPT*COR-XPT*SIR 

26230 

500 

CONTINUE 

26240 

RETURN 

26250 

END 

26260 

FUNCTION  A0N2PI (ANGIN) 

26270 

ANG=ANGIN 

26280 

IF(ANG)  1,2,3 

26290 

1 

ANG=ANG+360. 

26300 

IF(ANG) 1,2,2 

26310 

4 

ANG=ANG-360. 

26320 

3 

IF (ANG-360 . ) 2,4,4 

26330 

2 

A0N2PI=ANG 

26340 

RETURN 

26350 

END 

26360 

SUBROUTINE  LININT (Xl , Y1 ,X2 , Y2 , X , Y , ] 

26370 

ISTAT=-1 

26380 

DIST=0 . 0 

26390 

A1=Y2-Y1 

26400 

B1=X1-X2 

26410 

Cl=Xl*Al+Yl*Bl 

26420 

A2=DY 

26430 

B2=-DX 

26440 

C2=X*A2+Y*B2 

26450 

RAD=A2*Bl-Al*B2 

26460 

IF (RAD . EQ. 0 . ) RETURN 

26470 

XT=(B1*C2-B2*C1)/RAD 

26480 

YT=- (A1*C2-A2*C1) /RAD 

26490 

YMAX=AMAX1 (Yl,Y2) 
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26500 

26510 

26520 

26530 

26540 

10 

26550 

20 

26560 

30 

26570 

40 

26580 

26590 

26600 

26610 

60 

26620 

26630 

26640 

26650 

26660 

26670 

26680 

26690 

26700 

26710 

26720 

26730 

26740C 

*** 

26750 

875 

26760 

26770 

26780 

26790 

26800 

26810 

26820 

26830 

26840 

26850 

26860 

26870C 

*** 

26880 

876 

26890 

10 

26900 

26910 

26920 

26930 

26940 

26950C 

*** 

26960 

877 

26970 

26980 

26990 

YMIN=AMINl (Y1 ,Y2) 

XMAX=AMAX1 (XI ,X2) 

XMIN=AMIN1 (Xl , X2 ) 

IF(XT-XMIN+0. 00001)  60,10,10 
IF(XT-XMAX-0. 00001)  20,20,60 
IF (YT-YMIN+0. 00001)  60,3'0,30 
IF(YT-YMAX-0. 00001)  40,40,60 
DIST=SQRT ( (XT-X) **2+(YT-Y) **2) 

ISTAT=1 

DIFF=-C1+A1*X+B1*Y 

DIST=SIGN ( 1 . ,DIFF) *DIST 

CONTINUE 

RETURN 

END 

SUBROUTINE  CURRNT ( ALPHA , PSI , RADIUS , CTAN , CRAD , CROT) 
COMMON/CMLONG/TOWLEN , BOTLEN , DUM2 ( 4 ) ,NBARL 

COMMON/ CMROUT/ SDAT (11,10,2) ,CUR (11 , 10 , 3) ,SANG(11,10) ,NANG(10) ,NSEG 
COMMON/CMSTAT/DUM3 (17) , ISEG 
COMMON/CMDISP/DUM ( 5 ) , TLEN , TLEN2 , TLEN3 
COMMON/CMS I Z E/CGAFT , CGFWD , DUM1 ( 4 ) 

DIMENSIONXPOSN (3) ,CTI(3) ,CTO(3) ,CTA(3) ,CRA(3) 

XPOSN (1 ) =0 . 0 

XPOSN (2) =CGFWD 

XPOSN ( 3 ) =-CGAFT-BOTLEN 

WRITE (5,875) ALPHA, RADIUS , PSI 

FORMAT(/’  ALPHA' ,F7.3, ’ RADIUS ', F7 . 0 , ' PSI',F7.3) 

DO  10  N= 1 , 3 

CALL  BARDISl ALPHA, RADIUS, PSI, XPOSN (N) , Dl , D2 , PRAD , PTAN , 

&ISUBO, ISEGO) 

ISUBM=ISUBO-l 

CTI  (N) =CUR(ISUBO, ISEGO, 1 ) *PTAN+CUR ( ISUBM , ISEGO, 1)  * (1 .-PTAN) 

CTO (N) =CUR ( ISUBO , ISEGO ,2) *PTAN+CUR (ISUBM , ISEGO , 2) * ( 1 . -PTAN) 

CTA (N) =CTO (N) *PRAD+CTI (N) * ( 1 . -PRAD) 

CRA (N) =CUR (ISUBO , ISEGO , 3) *PTAN+CUR ( ISUBM , ISEGO , 3 ) * ( 1 . -PTAN) 

IF (ISEGO. EQ . ISEG) GOTO10 

CHSN=SIGN (1.0, SANG (1, ISEG) ) *SIGN ( 1 . 0 , SANG ( 1 , ISEGO) ) 

CRA (N) =CRA (N) *CHSN 

WRITE (5, 876) CTO (N) ,CTI (N) ,CTA(N) ,CRA(N) 

FORMAT ('  CTO  1 , F5 . 1 , ' CTI',F5.1,'  CTAN',F5.2,'  CRAN',F5.2) 

CONTINUE 

CTAN=CTA ( 1 ) 

CRAD=CRA(1) 

CLATB=CRA(2) *COS (PSI) +CTA (2) *SIN(PSI) 

CLATS=CRA ( 3 ) *COS (PSI) +CTA(3) *SIN (PSI) 

CROT= (CLATB-CLATS) /TLEN 
WRITE (5, 877) CTAN, CRAD, CROT 

FORMAT ( ' CURTAN’ ,F5.2, ’ CURRAD’ ,F5 . 2 , ' CURROT’ ,F5 . 2) 

RETURN 

END 

SUBROUTINE  WINDFO (ALPHA , PS I , FWLON , FWLAT , FWROT ) 
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27000 

27010 

27020 

27030 

27040 

27050 

27060 

27070 

27080 

27090 

27100 

27110 

27120 

27130 

27140 

27150 

27160 

27170 

27180 

27190 

27200 

27210 

27220 

27230 

27240 


COMMON/CMWIDE/TOWWID , BOTWID , NBARW 

COMMON /CMLONG/TOWLEN , BOTLEN , CGTOW , CGBOT , TOWK , BOTK , NB ARL 

COMMON/CMCHAR/TOWDRF , BOTDRF , TOWBC , BOTBC 

COMMON /CMSTAT/DUM3 (17) ,ISEG 

COMMON/CMCNST/RHO ,GRAV , PI 

COMMON/CMWIND/SWIND (10) ,DWIND (10) 

AX= (20 . -TOWDRF) *TOWWID 

AY= (15. -TOWDRF) * ( TOWLEN+BOTLEN) 

RHOA=RHO/800./2. 

ALREL=DWIND ( ISEG) *PI/180 . -ALPHA 
WSR=SWIND ( ISEG) *COS (ALREL) 

WST=SWIND (ISEG) *SIN (ALREL) 

WX=WSR*SIN (PSI) +WST*COS (PSI) 

WY=WSR*COS (PSI) -WST*SIN (PSI) 

CDX=1 . 0*RHOA*AX 
CDY=1.0*RHOA*AY 
CMX=0.0*RHOA*AX 
CMY=0.0*RHOA*AY 
WX=WX*ABS (WX) 

WY=WY*ABS (WY) 

FWLON=CDX*WX 

FWLAT=CDY*WY 

FWROT=CMX*WX+CMY*WY 

RETURN 

END 


